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ABSTRACT 

We present a novel pair of numerical models of the interaction history between the 
Large and Small Magellanic Clouds (LMC and SMC, respectively) and our Milky Way 
(MW) in light of recent high precision proper motions from the Hubble Space Tele- 
scope (jKallivavalil et alJl2QQ6alibh . These proper motions imply that the Magellanic 
Clouds (MCs) are moving ^80 km/s faster than previously considered. Given these 
velocities, cosmological simulations of hierarchical structure formati on favor a scenario 
where the MCs are currently on their first infall towards our Galaxy (Besla et alJl2007|; 
iBovlan-Kolchin eF aL 2011; Busha et al. 2011). We illustrate here that the observed 
irregular morphology and internal kinematics of the Magellanic System (in gas and 
stars) are naturally explained by interactions between the LMC and SMC, rather 
than gravitational interactions with the MW. These conclusions provide further sup- 
port that the MCs are completing their first infall to our system. In particular, we 
demonstrate that the Magellanic Stream, a band of HI gas trailing behind the Clouds 
150 degrees across the sky, can be accounted for by the action of LMC tides on the 
SMC before the system was accreted by the MW. We further demonstrate that the 
off-center, warped stellar bar of the LMC and its one-armed spiral, can be naturally 
explained by a recent direct collision with its lower mass companion, the SMC. Such 
structures are key mor phological characteristics of a cl ass of galaxies referred to as 
Magellanic Irregulars (jde Vaucouleurs fc Freeman! 119721) . the majority of which are 
not associated with massive spiral galaxies. We infer that dwarf-dwarf galaxy inter- 
actions are important drivers for the morphological evolution of Magellanic Irregulars 
and can dramatically affect the efficiency of baryon removal from dwarf galaxies via 
the formation of extended tidal bridges and tails. Such interactions are important 
not only for the evolution of dwarf galaxies but also have direct consequences for the 
buildup of baryons in our own MW, as LMC-mass systems are believed to be the 
dominant building blocks of MW-type halos. 

Key words: galaxies: interactions — galaxies: kinematics and dynamics — galaxies: 
evolution — galaxies: irregular — Magellanic Clouds 



1 INTRODUCTION 
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The Large Magellanic Cloud (LMC) is the prototype for 
a class of dwarf galaxies known as Magellanic Irregu- 
lars. Like the LMC, these galaxies are characterized by 
being gas rich, one-armed spirals with off-center bars 
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(|de Vaucouleurs Freemanlll972l ). Although there are nu- 
merous examples of Magellanic Irregulars in our Local Vol- 
ume, they are rarely found about massive spirals. This has 
been confirmed by recent studies of the frequency of LMC 
analogs about Milky Way (M^y) type gjalaxies in th e SDSS 
DR7 catalog (iLiu et al.ll2oTll: iTollerud et all 1201 ih. B ased 



19941: iGardiner Noguchil Il996l: bekki Chibal l2005l: 



on similar statistics, de Vaucouleurs F reeman (1972) con- 
cluded that the LMC has necessarily experienced little or 
no distortion due to interactions with the MW, and so its 
characteristic asymmetric features must owe to some other 
process. 

The idea that the LMC's evolution has not been dic- 
tated by interactions with the MW is given further cre- 
dence by the distance morphology relationship exhibited 
by MW and M31 satellites, whereby gas rich satellites 
are located at larger galactocentric radii than gas poor 
spheroidals. The Magellanic Clouds (MCs), at a mere 50- 
60 kpc away, are nota b le ex ceptions to this relationship, 
leading Ivan den Berghl (| 20061 ) to describe them as inter- 
lopers in our system. Along the same lines, recent stud- 
ies indicate that the LMC is much b luer in color relative 
to analogs in its m agnitude range (|Tollerud et al.1 l201ll : 
I James &: Ivorvllioilh . This fact is difficult to reconcile with 
the expected gas loss and quenching of star formation the 
LMC should have incurr ed if it were indeed a long-term 
companion of the MW (Grcevich &: Put man 2009 ). These 
conclusions are further supported by rec ent proper motion 
measurements (jKallivavalil et al.1 l2006al lbh. which indicate 
that the LMC is moving ^80 km/s faster than previously 
believed ([Gardiner Nog uchi 1996). Given the measured 
energetics of the LMC's orbit toda y, backward orbital inte- 
gration schemes (iBesla et al.ll2007|l and statistics froni large 
scale cosmol ogical simulations jBovlan-Kolchin et aI1l201ll : 
iBusha et al .l 2011) indicate that the LMC is likely on its ffist 
infall towards the MW. Consequently, the MW cannot have 
been the driver of its morphological evolution. 

In this study we ask the following: if not interactions 
with the MW, then what is the origin of the asymmetric 
appearance of the LMC and what is its connection to Mag- 
ellanic Irregulars in general? 

Notably, the LMC has a nearby companion, the Small 
Magellanic Cloud (SMC). I n fact, many M agellanic Irreg- 
ulars also have companions (Odewahn"l994'), although the 
frequency of such configurations is debated (Wilcots 2009^ 
Particularly striking exam ples include the Mage llanic Irreg- 
ular galaxies NGC 4027 (Fhookun et al.l Il992l ) and NGC 
3664 (jwilcots k Prescottll20ol )7 both have a low mass com- 
panion to which each is connected by a bridge of gas. The 
LMC and SMC are also connected by a bridg e of HI gas, 
known as the Magellanic Bridge (|Kerrl Il957f ). suggesting 
that interactions between dwarf pairs may hold clues to un- 
derstanding the current morphology of Magellanic Irregular 
galaxies. 

In addition to the Magellanic Bridge, the MCs 
are associated with both leading and trailing streams 
of gas, referred to as the Leading Arm and Magel- 
lanic Stream, respectively. The Magellanic S tream extends 
over 150 degrees across the so u thern skv ([Nidever et aD 
20101. l2008l: IPutman et all l2003l : IWannier WrixonI Il97 J 
"l974ll" and has been traditionally mod- 
eled as the product of MW tides (Murai & Fu iimoto 
Il980l : iHeller Rohi^ll994l : iLin et al.1 [1995. : .Gardiner et al 



Connors et al. 2005; Mastropietro et al. 2005; Ruzicka e t al 
2009, 2010) and ram pressure stripping (Mastropietro et al 
[2OO5I : iMoore Davislll994l ). A purely hydrodynamic solu- 
tion cannot pull material forward to explain the Leading 
Arm Feature, meaning that MW tides must be invoked in 
some form in all of these models. However, on a first in- 
fall, MW tides are negligible until very recently; it is thus 
difficult to reconcile the new proper motions and updated 
orbits with the formation of the Magellanic Stream, Bridge 
and Leading Arm in the context of t he existing scenarios. 

Alternativelv. iBesla et aD (|2010) (hereafter BIO) intro- 
duced a model to explain the observed large scale gas mor- 
phology of the Magellanic System thr ough tidal interaction s 
between the LMC and SMC (see also, IPiaz Bekkill2011al ). 
Because MW tides are not responsible for removing mate- 
rial from the system, this picture is consistent with a first 
infall scenario. In this model, the Magellanic Bridge, Arm 
and Stream are hypothes ized to be analogs of the classical 
iToomre k Toomre' (l972) tidal bridge and tail scenario and 
should be commonly found about interacting pairs/groups 
of dwarf galaxies. 

Here we explore whether interactions between the MCs 
can also account for the internal morphology and kinematics 
of the LMC and therefore shed light on the dynamical state 
of Magellanic Irregulars more generally. 

In particular, the nature of the LMC's off-centered stel- 
lar bar has been a long standing puzzle, as it is not present 
in any other tracer of the interstellar medium (ISM); it is 
neither apparent in the HI gas disk nor a site of active star 
formation as traced by Ka emission. Strong bars in more 
massive galaxies serve to funnel gas towards the center; 
streaming motions and characteristic "S-shaped" isoveloc- 
ity contours are thus evident in their gas velocity fields. 
While weak large scale streaming motions along the bar 
may be evident in the LMC HI velocity field (K im et al.l 
Il998h , the expected "S-shaped" isovelocity contours are not 
present Q Interestingly, this is also true of many other Magel- 
lanic Irregulars (Wilcots 2009): bars in these systems do not 
appear to strongly affect the underlying gas distribution. 
There is also evidence that the bar may be warped rela- 
tive to the LMC dis k plane (|Subramaniaml [2OO3I : iLah et ahl 
I2005I : Ik oerwer 2009 1) . Using r elativ e distance measurements 
to Cepheids, Nikola ev et al. liooi) 

concluded that the bar 
is in fact located ^0.5 kpc in front of the main disk . For this 
reason it has been described as a "levitating" bar. IZaritskvl 
(2004) suggests that this may be a result of viewing a tri- 
axial stellar bulge that is embedded in a highly obscuring 
thick disk. Along the same lines, Zhao & Evans (2000) pos- 
tulate that the off-centered bar is an unvirialized structure, 
inclined relative to the plane of the LMC disk by as much 
as 25 degrees in order to explain the microlensing optical 
depth observed towards the LMC. Clearly, the nature of the 
LMC's bar is an ongoing subject of debate. 

We posit here that a recent direct collision between 
the LMC and SMC has left the LMC with a warped, off- 



^ iKim et al] (|l998l ) comment on the existence of a distorted S- 
shap ed isovelocity con tour across the LMC's minor axis. How- 
ever, lOlsen MassevI (|2007h find that this feature straightens 
out when the higher proper motions are accounted for. 
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centered stellar bar and pronounced one-armed spiral. We 
further claim that such asymmetric structures are character- 
istic of Magellanic type galaxies undergoing minor mergers. 
In this study we illustrate that such a scenario is consistent 
with a first infall towards our MW and can simultaneously 
explain both the morphology and kinematics of the LMC 
as well as the large scale gas morphology of the Magellanic 
System. Thus, Magellanic Irregulars with nearby compan- 
ions should also be associated with faint extended gaseous 
tails and bridges. As in the Magellanic System, such features 
could hold ^50% of the baryonic mass of the original system, 
indicating that dwarf-dwarf tidal interactions are an impor- 
tant mech anism for the loss of b aryons in low mass systems 
(see also, IP'Onghia eraP l2009l ). as a con sequence of res- 
onan t interactions between spinning disks ( D'Onghia et al.l 
I2OIQI ). 

We stress that the goal of our study here is not to repro- 
duce every detail of the Magellanic System, as we have not 
conducted a complete parameter search of all the possible 
orbital configurations, mass ratios and gas fractions, which 
influence the final outcome. The aim of this investigation is 
rather to determine which of the observed peculiarities of 
the Magellanic System can be directly linked to interactions 
between the MCs. 

Moreover, our work has broader implications for under- 
standing the properties of accreted satellites. Minor merg- 
ers are frequent events that shape galaxies and their ha- 
los; however, little attention has been given to the accre- 
tion of binary pairs or groups of smaller galaxies (but see, 
'D'Onghia Lakel l2Q08: Sales et al. 2007). This study repre- 
sents a first step towards understanding the morphological 
evolution and gas loss rates of such galaxies immediately 
after their capture by a massive host. LMC mass objects 
are expec ted to be the p rimary building blocks of MW type 
galaxies (|Stewart et al.| [2008). making this study of direct 
relevance to our understanding of the evolution of the MW. 

In this paper we begin by outlining our methodology 
and introducing two possible models for the interaction his- 
tory of the MCs, one of which invokes a recent direct collision 
between the MCs. In the subsequent sections we discuss the 
resulting large scale gas structure and internal structure and 
kinematics of the LMC. The results for the SMC and the ex- 
pected stellar counterpart to the Magellanic Stream will be 
presented in future work. 



2 METHODOLOGY 

We follow the general method outlined in BIO to set up the 
initial galaxy models and orbits in order to reproduce the 
observed large scale gaseous structure of the Magellanic Sys- 
tem. Details about our numerical methods, initial conditions 
and chosen orbital parameters are described below. 



2.1 Numerical Methods 

All of the numerical simulations performed in this work 
use the N-bod y smoothed-pa rticle hydrodynamics (SPH) 
code. Gadgets (|SDringelll2005l V The Gadgets code incorpo- 
rates a subresolution mul tiphase model of the ISM th at in- 
cludes radiative cooling (Springel & Hernquist 200S), and 
incorporates a fully conservative approach to integrating 



the equations of motion (jSpringel &: HernquistI 1200^ 1. Star 
formation from the cold phase (i.e. all cold gas - no dis- 
tinction is made between atomic and molecular compo- 
nents) follows a Schmidt volume density law psFR oc 
(with N — 1.5) that is normalized to approximate the 
star formation rate of the MW. A local threshold volume 
density cutoff of 0.004 Mopc"^ is adopted, below which 
stars do not form. As pointed out by many authors (e.g.. 



Kuhlen et aLll201ll:lHopkins et al.feoilUGnedin Kravtsovl 



20 id : iRobertson Kravtsovl l2008l ) , such a star formation 
prescription is likely inappropriate for dwarf galaxies. We 
discuss the implications of our adopted prescriptions to our 
results in § 16.11 

Stellar feedback in the form of galactic winds is not 
employed in our simulations; however, we comment on the 
relative importance of outflows to the formation of the Mag- 
ellanic Stream in Appendix [B] 

We note that the reliability of SPH for cosmologi- 
cal simulations ha s recently been called into question by 
Vogelsberger et al." ('2011'); 'Siiacki et al." feOllV "Ke res et al.l 
(2011); Torrev et al. (2011); Bauer & Springel (20li). How- 
ever, comparisons between SPH and calculatio ns done with 
the moving mesh code Arepo (|Springell l2010h , show good 
agreement for applications involving galaxy collisions, at 
least when the subresolution model mentioned above is used 
to represent star-form ing gas (Hayward e t al. 2011, in prep). 
The tests done by e.g. lSiiacki et al.l(|201lh indicate that SPH 
can fail when applied to situations in which gas in very dif- 
ferent phases are in motion relative to one another. The use 
of an effective equation of state to describe the ISM effec- 
tively circumvents this issue because the different phases of 
the gas are not modeled explicitly. 

2.2 Initial Conditions 

The initial conditions for the construction of the LMC and 
SMC galaxies used for all models are outlined in Table [1] As 
in BIO, the total initial mass of the LMC is determined using 
current halo occupation models to relate the observed stellar 
mass of the L MC to its origina l halo mass before infall into 
the MW halo (I Cuo et al1l2010h . Reflecting their stellar mass 
ratio, the SMC is then chosen to be 10 times less massive 
than the LMC. Consequently, the MCs are modeled here to 
have infall masses an order of magnitude larger than em- 
ployed in previous models. The number of particles of each 
component (gas, stars, dark matter) are chosen such that 
the mass resolution per particle of a given type is roughly 
the same in both galaxies. 

The SMC is modeled with an extended gaseous disk 
with a scale length S times that of the stellar component. 
Much larger ratios are common for isolated dwarfs found 
in voids (Kreckel et al. 2011), and neutral hydrogen obser- 
vations of SMC-like dwarfs with t he Westerbork Synthesis 
Radio Telescope bv ISwaters et al. I (12002.) find HI disk scale 



lengths ranging from 1.4-4.5 kpc ( Connors et al.ll2005[ ). Our 

adopted scale length of S.S kpc is consistent with the upper 
end of the observed range. 

The LMC is modeled with gas and stellar disks with 
the same scale length, rather than with an extended gaseous 
disk. In reality, the interaction with the ambient hot gaseous 
halo of the MW would serve to truncate the LMC's extended 
gas disk. The scale height of the stellar disk is taken as 0.2 of 
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the disk scale length. The modeled scale height of the LMC's 
stellar disk is thus initially zo= 0.34 kpc (the observed value 
today is Rdisk = 1.4 kpc and zo= 0.27 kpc; van der Marel 
et al. 2002). The gaseous disk height is determined by self- 
gravity and the pressurization of the ISM, as prescribed b y 
the chosen effective equation of state (jSpringel et al.ll2005l l. 

The dark matte r halos of the LMC and SMC follow 
Hernquist potentials (jHernquist 1990 ). The scale radius for 
the Hernquist potential (th) is related to the scale radius of 
the c orresponding NF W halo (rs — R200 1 C ) ([Navarro et al.l 
I1997I ) as described in lSpringel et all l|2005h : 



r„ = rs^2[\n{l + C)-j-^j, (1) 

where C is the concentration parameter. Values for C, rs 
and rH are listed in Table [H 

The MW is modeled as a static NFW potential with a 
total mass of 1.5 x 10^^ M©, C = 12, virial radius of Rvir = 
300 kpc, and R200 = 220 kpc (radius where the average den- 
sity is 200 times the critical density of the Universe). As in 
BIO, dynamical friction from the MW halo is not explicitly 
accounted for, but is expected to have little impact on the 
orbit in a first passage (see Besla et al. 2007, Figure 4). Dy- 
namical friction between the MCs, on the other hand, plays 
a much more important role in their orbital evolution and is 
captured explicitly by modeling these two galaxies with live 
dark matter halos. 

The resulting rotation curves for the MCs are plotted in 
Figure [U The initial SMC rotation curve peaks at Vrot = 60 
km/s at 3 kpc from the center, as expected from HI kinemat- 
ics (Stanimirovic et al. 2004); the SMC is initially a well- 
behaved disk galaxy. The initial simulated LMC rotation 
curve peaks at Vrot — 95 km/s, which is within the observed 
range (IStavelev-Smith et al.1 l2003al : lOlsen MassevI l2007l : 
Ivan der Marel et al.ll2002l ). 



2.3 Orbit of the SMC about the LMC and 
Definition of Models 

Following the method outlined in BIO, the MCs are evolved 
as an isolated interacting binary pair over a period of ^7 
Gyr since the SMC first crossed within i?2oo = 117 kpc of 
the LMC. 

The SMC is placed on an eccentric orbit about the LMC 
(ecc = 0.7). Higher orbital eccentricities for the SMC result 
in fly- by encounters between the MCs, while lower values 
cause the SMC's orbit to decay too quickly. 

The simulation is stopped at characteristic points in 
time, defining two models for the orbital history of the SMC 
about the LMC, referred to as Model 1 and Model 2. Model 

1 is stopped after 5.1 Gyr and Model 2 after 5.9 Gyr. Thus, 
Model 1 and 2 differ based on the number of passages the 
SMC has completed about the LMC. In Model 1, the SMC 
has completed 2 passages about the LMC, whereas in Model 

2 it has completed 3. The stopping times are chosen such 
that 1 Gyr after this time, the LMC will have travelled from 
a distance of 220 kpc (i?2oo for the MW) to its current loca- 
tion and the SMC will have completed the desired number 
of orbits about the LMC. 

The choice of these two models is motivated by the 



overarching goal of this study to assess the role of interac- 
tions between the MCs to their evolution. Given the chosen 
L/SMC mass ratio of 1:10, it is unlikely that the SMC could 
have survived more than 3 passages about the LMC, mak- 
ing Model 2 a maximal interaction scenario. Model 1 is very 
similar to the solution presented in BIO - the analysis of 
such a model is a direct extension of the BIO work. Tidal 
forces between the MCs have been acknowledged as playing 
an important role in the formation of the Magellanic Stream 
in many previous studies. In fact, in IConnors et aD (|2006l ). 
the tidal force from the LMC on the SMC dominates over 
MW tides for most of the SMC's orbit. But to explain the 
Leading Arm Feature and extent of the Magellanic Stream, 
MW tides have been invoked in all of these studies. Instead, 
here and in BIO we illustrate how such extended structures 
can form without requiring the MCs to complete an orbit 
about the MW. 

The orbit of the SMC about the LMC in Model 1 and 
Model 2 is plotted in the top panel of Figure O The black 
line indicates the evolution of the system in isolation (no 
MW potential) and is continued 1 Gyr past the respective 
stopping point for each model. The red line shows how the 
orbit of the SMC is modified if instead the binary pair is 
captured by the MW after the stopping point. 

In Model 1, MW tides work to keep the LMC and SMC 
further apart than they would have been in isolation. In 
Model 2 the opposite occurs; the MW's gravitational pull 
forces the SMC to colhde directly with the LMC. The dis- 
tinct outcomes occur because of differences in the SMC's 
separation from the MW relative to its separation from the 
LMC at the pericenter of its orbit about the MW. At peri- 
center, the SMC is closer to the LMC in Model 2 than in 
Model 1 (see Figure (H). 

The bottom panels of Figure [2] show the gas distribu- 
tion of the L/SMC at the stopping points for each model as 
contours plotted over the stellar distribution. After 5.1 Gyr 
(stopping point for Model 1) the SMC is at the apocenter 
of its second orbit about the LMC. After 5.9 Gyr (stopping 
point for Model 2) the SMC is just completing its second 
orbit. As outlined in BIO, a tidal bridge and extended tail 
forms as a result of the action of LMC tides on the SMC; 
these features form before the system is captured by the 
MW. 

2.4 Orbit of the MCs about the MW 

At the stopping time (5.1 Gyr for Model 1 and 5.9 Gyr for 
Model 2), the isolated MC pair is placed at i?2oo = 220 kpc 
from MW's galactic center, as illustrated in Figure [3] (left 
panels). i?2oo is chosen as the starting radius because the 
MW tidal field does not distort the orbit of the SMC rel- 
ative to the isolated orbit until well within that radius; it 
takes 500 Myr for the red and black lines in the top pan- 
els of Figure [2] to deviate after the stopping point (i.e. 500 
Myr after they cross i?2oo)- As such, the overall interaction 
history of the MCs is well-described by the isolated system 
before this point. 

The galaxies travel to their current locations on orbits 
consistent with the HST proper motions for the LMC, as 
indicated in Figure [3] (right panels). This takes 1 Gyr in 
both models since the LMC's orbit about the MW is roughly 
the same in both cases - it is the SMC's orbit that differs. 
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Figure 2. Top Panel: The separation of the SMC from the LMC is plotted as a function of time for Model 1 (left) and Model 2 (right). 
The black line denotes the orbital history when the LMC and SMC are modeled as an isolated binary pair: the MW's potential is not 
included. The red line indicates how the separation between the MCs is modified after the pair first cross within R200 = 220 kpc of the 
MW. This occurs after 5.1 Gyr for Model 1 and after 5.9 Gyr for Model 2; these times are referred to as the stopping times for each 
model. Bottom Panel: The gas column density of the isolated LMC-SMC system is projected in the binary orbital plane as contours over 
the stellar distribution at the stopping time for Model 1 (left) and Model 2 (right). Gas contours span a range of lO-*^^ — lO^*-* cm~^, 
where each contour represents an increase in column density by a factor of 1.5. These images thus depict the gas distribution of the 
simulated system before it is influenced by the MW's gravitational potential. An extended tail of gas is stripped from the SMC and a 
bridge of gas connects it to the LMC in both models (see also Figure 1 of BIO). 
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Table 1. L/SMC Initial Conditions 



Property 


LMC 


SMC 


M* (Mo) 


2.5 X 10^ 


2.6 X 10^^ 


Mgas (Mo) 


1.1 X 10^ 


7.9 X 10^ 


^gas 


0.3 


0.75 


Mtotal (Mo) " 


1.8 X 10^1 


2.1 X 10^° 


R200 (kpc)^ 


117.1 


57.1 


C 


9 


15 


(kpc)^ 


13.0 


3.8 


th (kpc)^ 


21.4 


7.3 


Stellar Disk scale length (kpc) 


1.7 


1.1^ 


Gas Disk scale length (kpc) 


1.7 


3.3 


Gravitational Softening Gas/Stars (kpc) 


0.1 


0.1 


Gravitational Softening Halo (kpc) 


0.29 


0.29 


Nstars 


10^ 


10^ 


Ngas 


3 X 10^ 


3 X 10^ 


Nhalo 


10^ 


10^ 




0.3 


0.3 



"Note that the initial stellar mass and stellar disk scale length chosen for the SMC deviate from the values adopted in BIO. Here the 
disk is chosen to be more extended in order to increase the number of stars removed by LMC tides. Other changes in parameter values 

are minor and reflect attempts to match various observed mass constraints for the MCs (see Table [2]). 
^The gas fraction relative to the total disk mass (stars + gas). The ^as fract ions of isolated dwarf galaxies are known to be large, e.g. 

iGeha et al.l (|2006h 

^The total mass of the LMC /SMC at infall is determined using the observed stellar ma ss of the L MC(SMC) M* = 3 x 10^ Mo 
(3 X 10^ Mo) (|van der Marel et al.ir2002l : IStanimirovic et al.ll2004h and the relations from lCuo et al.l ('iolO). The total halo mass used 
to deflne the Hernquist dark matter proflle is then M^alo = ^total - - Mgas- 
^The radius where the average enclosed density is 200 times the critical density of the universe 
^The scale radius for the NFW proflle (rs), which is used to deflne the scale radius of the Hernquist proflle (th), following 

Springel et al. (2005). 

•^The effective equation of state parameter, q, deflnes the pressurization of the ISM following (jSpringel et al.lEoOSl V 



The Galactocentric position and velocities of the MCs 
are plotted in Figure |4] as a function of time since they first 
crossed within R200 of the MW. Also plotted are the relative 
positions and velocities between the MCs (orange line). 

In Model 2 the SMC completes an additional passage 
about the LMC since entering the virial radius, versus in 
Model 1. This additional passage results in a direct collision 
between the MCs and the formation of a new bridge. Tidal 
bridges and tails are form ed at each pericentr ic passage of 
the SMC about the LMC (jToomre Toomrel [T972V Thus, 
in Model 2, the bridge connecting the MCs will have formed 
^100 Myr ago, during this direct collision (separation ap- 
proaching zero). 

2.5 Comparison of Modeled Orbits with Data 

For each model, the initial velocities and positions of the 
LMC and SMC at R200 = 220 kpc from the MW and 
their final values today are summari zed and compared to 
data from iKalliyavalil et alJ (|2006ah (hereafter Kl) and 
iKallivavalil et al J (|2006bD (hereafter K2) in Table El The 
proper motion error space for the Kl HST proper mo- 
tion measurements of the LMC is indicated in Figure \5\ 
Over-plotted are various other measurements for the LMC's 
proper motion and the simulated Model 1 and 2 results. The 
final LMC velocities and positions are designed to be within 
la of the observations in both models. The differences be- 



tween the two models are the orbital parameters for the 
SMC. 

In Model 1 the SMC velocity is significantly larger than 
that indicated by the HST proper motions of K2. In Model 
2, the SMC velocities are in better agreement, however the 
separation between the LMC and SMC is smaller than ob- 
served (by about 10 kpc). As such, the line-of-sight velocity 
and proper motions for the SMC are also different than ob- 
served. 

While the SMC velocities and positions are not perfect 
matches to the observations, it is unlikely that significantly 
new insight would be gained as to the physical processes 
at work if an exact solution were found. Slightly different 
choices of orbital parameters and timing in the orbit can 
change the SMC's final position and velocity, but not the 
physical picture. This is practically illustrated by compar- 
ing the resulting large-scale gas distribution in Model 1 and 
2 (see § [3]); despite differences in the SMC orbital proper- 
ties, the same overall scenario has produced similar global 
features (i.e. a Leading Arm, Bridge and Stream). To match 
the exact properties of the Magellanic System, a more de- 
tailed study, varying orbital parameters, L/SMC mass ra- 
tios, MW mass, etc, is required; this is beyond the scope of 
the present study. 

We have obtained another epoch of data with WFC3, 
resulting in an average time baseline of 7 years (Kallivayalil 
et al. in prep). These new data are expected to reduce the 
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Figure 3. Projected gas column densities in the YZ galactocentric plane for the simulated Magellanic System when the LMC first crosses 
within i?200 of the MW (time=0, left) and today (time~1.0 Gyr, right). The circle indicates the location of i?200 = 220 kpc. The results 
for Model 1 are plotted in the top row and Model 2 is on the bottom. The LMC's(SMC's) orbital path is denoted by the solid(dashed) 
line. 



errors on the proper motions by a factor of 3, potentially 
narrowing parameter space the SMC's error space. We note 
that, within this error space, the exact choice of LMC and 
SMC velocity today will not alter the physical picture pre- 
sented in this work, which is that tidal interactions between 
the two Clouds are the main driver for their morphological 
and kinematic evolution. 



3 LARGE SCALE GAS MORPHOLOGY 

The resulting large scale gas distributions in Models 1 and 
2 are shown in a Hammer- Aitoff projection in Figure [6l In 
both models, the final gas distribution can be described as 
an extended tail, a leading component and a bridge of gas 
connecting the two galaxies. As such, the main components 
of the Magellanic System are reproduced by both models. 
Moreover, in both cases the simulated stream stretches ^150 
degrees across the sky, as observed (Nidever et al.ti2QQ8. ). 



In Figure [71 the simulated stream is plotted in Mag- 
ellanic Coordinates, a variation of th e galactic coordinat e 
system where the Stream is straight ([Nidever et al-llioOSl ). 
In both models, the simulated stream deviates away from 
the projected location of the past orbits on the plane of the 
sky, as expected ac cording to the r ecent proper motions (see 
e.g.. Figure 8 in lBesla et al.1 120071 ). The deviation is a nat- 
ural result of the proposed stream formation mechanism. It 
occurs largely because the Stream is removed in the binary 
LMC-SMC orbital plane by LMC tides. This binary plane is 
not parallel to the LMC-SMC-MW orbital plane, thus the 
Stream is not coincident with the orbit of the MCs about the 
MW. A second factor is the orientation of the SMC's disk; 
the location of the simulated stream can be modified by 
changing this angle. In both of these models the SMC disk 
is initially oriented 90 degrees with respect to the SMC- 
LMC orbital plane. The deviation between the simulated 
stream and the orbits is more pronounced in Model 1 than 
in Model 2. However, this could be altered if the SMC disk 
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Figure 4. The Galactocentric position (top) and velocities (bottom) of the LMC (green) and SMC (black) are plotted as a function of 
time since the MCs first crossed within R200 = 220 kpc of the MW. The relative separation and velocity between the MCs is plotted in 
orange. The results of Mod el 1 are plot ted in the left column and those of Model 2 are in the right column. The Galactocentric velocities 
determined by (|Kallivavalil et al.ll20Q6al lU hereafter Kl and K2) are 378 ± 18 km/s for the LMC and 302 ± 52 km/s for the SMC (errors 
quoted are 1 a). In Model 1, the velocity of the SMC is higher than observed. In Model 2, the galactocentric position of the SMC is 
too small (53 kpc vs. 60 kpc); consequently the relative separation between the MCs is also too small (11 kpc vs. 23 kpc). Note that in 
Model 2, the separation between the LMC and SMC approaches zero ~ 100 Myr ago, indicative of a direct collision. 



Table 2. Initial and Final Orbital Parameters: Model 1 and 2 







Model 1 


Model 2 


Observed Today 


Galaxy 


Parameter 


At R200 


Today 


At R200 


Today 


Kl and K2 


LMC 


(x,y,z) (kpc) 


(35, 203, -63) 


(-1, -40, -25 ) 


(48, 198,-85) 


(-1, -42, -26) 


(-0.8, -41.5, -26.9) 




(vx,vy,vz) (kpc) 


(-14, -157, -29) 


(-72, -267, 250) 


(-17, -160, -29) 


(-82, -263, 249) 


(-87 ± 12, -268 ± 11, 252 ± 16) 




Vios (km/s) 




262 




259 


262 ± 3.4 




PM (W, N) (mas/yr) 




(-2.02, 0.44) 




(-2.03, 0.43) 


(-2.03 ± 0.08, 0.44 ± 0.05) 


SMC^ 


position (x,y,z) (kpc) 


(5, 243, -62) 


(18, -46, -46) 


(56, 193, -90) 


(6, -39, -35) 


(15.3, -36.9, -43.3) 




velocity (vx,vy,vz) (kpc) 


(6, -146, -70) 


(-88, -384, 246) 


(-51, -289, 88) 


(-66, -258, 198) 


(-87 ±48, -247 ±42, 149 ± 37) 




Vios (km/s) 




215 




201 


146 ± 0.6 



^ The SMC proper motions are not included in this table: since the line of sight velocities are not well-matched to the observations the proper motions 
can't be meaningfully compared to the data. 

were oriented differently initially and is not a physical dis- The structure of the Leading Arm Feature (LAF) is dis- 
tinction between the models; the magnitude of the offset is a tinct in each model. In Model 1 the LAF represents material 
tunable parameter. Note that this offset is not expected in a that was stripped from the SMC on earlier passages and cap- 
ram pressure solution for the Stream, as the material should tured by the LMC. Since this material is bound to the LMC 
be removed along the direction of motion (see Appendix[A|. it does not extend further than 50 degrees (see Figure [T]). 



© 2012 RAS, MNRAS 000, [HS] 



The Magellanic System 9 





1 1 1 1 1 1 1 1 1 1 




1 ' ' 




- □ P08 


X 


Kl (m^M^> 




; :4^vdM02 




vio ; 


0.6 






Model 1 - 






^^^-^ ■ 
^--.^ '\ 


I^'Odel 2 ; 

\ 


^0.5 






\ \ 


^0.4 


{ { { { ^ 

1 I I \^ 


b ] 


\ 1 

/ 1 
/ / 
/ / — 








/ / 
/ / 










0.3 








0.2 


- LMC 

7 1 1 11 1 1 1 1 1 




1 1 1 r 



Mw (mas yr-i) 



Figure 5. Concentric circles indicate the 4cr error space for the Kl proper motion error space for the LMC (where the mean value is 
indicated by the X). The asterisk indicates the mean of all proper m otion estimates for the LMC prior to 2002 (^van der Marel et alJ 
I2QO2I V The circled dot indicate s the recent proper motion estimate bv lVieira et alJ (|2010l ) and the open square shows the reanalysis of 
the Kl proper motion data bv lPiatek et alJ (|2008l V The red triangle and blue square indicate the proper motion of the LMC today in 
Model 1 and 2, respectively. These values were chosen to closely match the Kl data. 



In Model 2, the LAF is better described as a tidal tail or 
loop, resulting from tidal stripping of the SMC on its most 
recent orbit about the LMC. The tidal tail gains energy and 
deviates to large r distances away from the SMC's orbit (e.g., 
IChoi et aLlboOTl l. The resulting simulated LAF in Model 2 
spans 80-90 degrees across the sky, which is larger than ob- 
served. Unlike the Stream, this material is leading to the 
MCs and so will experience a significant ram pressure head- 
wind. Consequently, its final appearance, position and an- 
gular extent on the plane of the sky cannot be well captured 
without including hydrod ynamic effects (see Appendix [Al 
iBinnev &: Fraternalill201lh . Model 2 does illust rate, however, 
that the observed 70 degree span of the LAF (|Nidever et al.l 
I2OIOI ) can be reproduced without invoking a previous orbit 
about the MW. 

The Magellanic Stream is observed to have a pro- 
nounced HI coluran density gradient along its exten t 
(jNidever et al.ll2010l : IPutman et al]l2003l : iBriins et alll2005l l . 
The maximum gas column density along the simulated Mag- 
ellanic System is determined from Figure [71 and plotted as a 
function of Magellanic Longitude in Figure [51 Both models 
underestimate the observed values, which are indicated by 
the solid red line (Briins et al. 2005). There are a number 
of possible explanations for this discrepancy. This problem 
could be alleviated if the gas reservoir in the SMC were 
depleted less efficiently at earlier times; for example, if the 
SMC's gas disk were initially less extended or if star forma- 
tion was not quite so efficient. Ram pressure stripping has 
falso not been modeled and could also increase the amount of 
gas removed as the MCs get closer to the MW. Furthermore, 
hydrodynamic instabilities are not we ll modeled with SPH 
(jAgertz et al.ll2007l : Isiiacki et al.ll201ll l. and so clumping of 
the gas is not captured in these simulations. This is a process 



that will clearly influence the resulting gas column densit^n 
the Stream (Bland-Hawthorn et al. 2007 : Nigra et al. 201QI 1. 

At the same time, there are notable consi s tencie s be- 
tween the models and the data. iNidever et al ] (|2010") find 
that the column density along the LAF is fairly constant 
along its ^-^70 degree span (^-^ 4 x 10"^^ cm~"^); this is true 
of Model 2. The column density in the bridge in Model 2 
also matches the observations: the maxi mum colurnn den - 
sity in the Bridge is 1.64 x 10^^ cm"^ (Briin s et ai]l2005l ): 
the simulated bridge column density for Model 1 is too low. 
Also in Model 2 the simulated column density of the SMC 
is hig her than that of the LMC, as observed ([Briins et al.l 
I2OO5I . the maximum column density is 5.45 x 10^^ cm~^ and 
9.98 X 10^^ cm"^ for the LMC and SMC, respectively). 

There is a well-d efined velocity grad i ent along the 
lengt h of the Stream (I Putman et al ] I2OO3I : INidever et al.l 
2008), ranging from 200 km/s near the Clouds to -400 km/s. 
The simulated results are plotted in Figure [H for both mod- 
els. The observed line-of-sight velocities along the MS (yel- 
low line) and the rest of the system are well-traced by Model 
1. Model 2 also reproduces the observed range of velocities, 
but the slope of the velocity gradient along the MS is not 
well-matched to the data. 

Given that the only difference between Model 1 and 
Model 2 is the SMC's orbital parameters, it is doubtful 
that this discrepancy in predicted velocities owes to miss- 
ing physics. Rather, a detailed search of the SMC's orbital 
parameter space will likely yield better matches for Model 2. 
It is possible that gas drag effects (not modeled here) may 
also modify the velocity profile - particularly in the LAF, 
where the velocities are currently too high in Model 2. 

No direct distance measures exist for the Magellanic 
Strea m, as no stellar c ounterpa rt has yet be en identified 
(fGuhathakurta L ReitzeL. 1998. 1 iJin Lvnden-Bell (2008 ) 
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Figure 6. Hammer-Aitoff projection of the total gas distribution of the simulated Mage llanic System (r ed scale) for Model 1 (top) and 
Model 2 (bottom) is plotted over an image of the MW (blue, white and brown colors; iMellingeJ l2009l ) . The orbital trajectory of the 
LMC(SMC) is indicated by the solid(dotted) white line. Various components of the Magellanic System are labelled, where LAF stands 
for Leading Arm Feature. 



present a geometrical method to determine distances along 
streams with well defined velocity gradients. Using this 
method t hey find the ti p of t he 100 degree long Stream 
defined in Putman et al. (2003) (i.e. not including the ex- 
tension recently described by Nidever et al. 2010) to be lo- 
cated at 75 kpc from the Galactic center. The line-of-sight 
distances of the gas in the simulated Magellanic System are 
plotted in Figure [101 The stream produced by Model 1 is 
generally closer (80-150 kpc) than that of Model 2 (80-230 
kpc). Both si mulated streams ar e furth er away than pre- 
dicted by the Ijin k Lvnden-Belll (|2008l ) method; however, 
gas drag and changes in the model parameters (such as in- 
creasing the MW mass) can alter the distance to the simu- 
lated stream. 



4 LMC MORPHOLOGY 

In this section we study in detail the resulting structure 
of the simulated LMC stellar and gaseous disks in our two 
models of the large scale gas distribution of the Magellanic 
System. 

Figure [TTI shows the LMC's stellar disk in Model 1 (left) 
and Model 2 (right) in our line-of-sight view. The RA and 
DEC coordinate grid is overplotted in green across the face 
of the disk. In both models the LMC disk is inclined ^35 
degrees with respect to the plane of the sky, as observed 
(i.e. despite the recent collision of the SMC in Model 2, the 
inclination of the LMC's disk remains unchanged). 

The Model 1 disk is fairly uniform and symmetric. In 
Model 2, however, there are perturbations induced in the 
LMC's stellar disk by the recent encounter with the SMC. 
In particular, there are significant distortions in the North- 
East. Only LMC stellar particles are plotted in these images. 
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Figure 7. The total gas distribution of the simulated Magellanic System for Model 1 (top) and Model 2 (bottom) is plotted in Magellanic 
Coordinates. The orbital trajectory of the LMC(SMC) is indicated by the solid(dashed) yellow line. The actual location of the Magellanic 
Stream is roughly traced by the solid white line. The white circle represents the observed radius of the LMC. 
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Figure 8. The maximum total gas column density of the simulated Magellanic System is plotted as a function of Magellanic Longitude 
for Model 1 (black line) and Model 2 (green thick line). The red lines indicate the observed maximum HI column density from the 
iBriins et al ] (l2005l ) data set for each marked region: MS stands for Mage llanic Stream and LA F for Leading Arm Feat ure. The various 
roman numerals refer to specific sections of the MS/LAF as defined in iPutman et al.l (|2003l ) and iBriins et al.l (|2005l ). Neither model 
reproduces the observed HI column density gradient in the MS. 



and so these structures are in the plane of the LMC disk and 
do not represent tidal debris from the SMC. Such structures 
are observed in deep observations of the perip hery of the 
LMC's disk (|de Vaucouleurs k FreemanI [l972l . Martinez- 
Delgado et al. in prep). 

The initial LMC disk was bar unstable, and so the bar 
feature in both models was present since the beginning of 
the simulation - it was not induced by external tidal per- 
turbations from the SMC or MW. Interestingly, in Model 2, 



the bar of the LMC is now off-centered relative to the disk, 
as observed. No such perturbations are observed in Model 
1: without a close encounter the LMC looks like a symmet- 
ric spiral disk galaxy and it is doubtful that such a galaxy 
would be classified as a Magellanic Irregular galaxy. 

In Figures [12] and [13] we take a closer look at the LMC's 
gas and stellar disk by deprojecting the disk from the line- 
of-sight frame into a Cartesian coordinate system centered 
on the LMC disk plane for both Model 1 and Model 2. Only 
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Figure 9. The line-of-sight velocities for the total gas distribution of the simulated Magellanic System are plotted a s a function of 
Magellanic Longitude for Model 1 (top) and Model 2 (bottom). The yellow line is a fit to the data of Nide ver et al The modeled 

line-of-sight velocities along the past orbit of the LMC(SMC) are plotted as the solid (dotted) white line. 



particles associated with the LMC are plotted. The images 
are centered on the peak of the stellar density distribution 
(i.e. the photometric center). 

In Model 1, the bar of the LMC is clearly visible in 
both the stellar and gaseous disks. This i s in d isagreement 
with the observed HI maps of iKim et aP (|l998l l: there IS no 
distinguishable bar structure in the observed global HI emis- 
sion that is comparable to the optical bar. The results for 
Model 2 illustrate the consequences of a recent (100-300 Myr 
ago) direct collision between the LMC and SMC. In this sce- 
nario, Model 1 represents the state of the LMC disk before 
the collision occurred, where the LMC is a symmetric barred 
spiral. After the collision, the bar has become off-centered 
with respect to the underlying disk and it has almost disap- 
peared from the gas disk of the LMC. The reason becomes 
clear when the disk is viewed edge-on along the x-axis. The 
bar has become warped by ^10-15 degrees relative to the 
LMC disk plane and is therefore inefficient at funneling gas 
in the way it could in Model 1. This warp o f the bar is less 
extrem e than that required by models such as I Zhao EvansI 
(|200Q| ) for the microlensing optical depth and may be consis- 
tent w ith the observations of S ubraman iam & Subramanian 
(|2009l ). The simulated offset bar is also consistent with the 
structure of the observed bar, which is described as a cigar 
shaped structure with dimensions of 1 x 3 kpc. 

From the edge-on view of the Model 2 gas disk, it is 
clear that LMC gas particles have been pulled out of the 



disk by the passage of the SMC through the LMC. This 
causes the appearance of a gaseous "arm" in the face-on 
view of the disk. Such "arms" are seen in the LMC gaseous 
disk (Kim et al. 1998) and are believed to be related to the 
Magellanic Bridge and Leading Arm Feature jNidever et al.l 
.2008 ). In our interpretation, at least one of these "arms" 
is extra-planar and located behind the LMC disk. The fact 
that LMC gas is removed from the disk towards the Bridge 
indicates that the formation of the Magellanic Bridge has 
been aided by hydrodynamic gas drag. It is not purely a 
tidal feature. The Bridge is known to be quite metal poor 
along two sightlines to wards early type stars (see Table IBll 
and Lehne r et al.l l2008Y However, a full census of the metal- 
licity across the Bridge does not yet exist. Regardless, it is 
clear that LMC gas cannot have contaminated the entirety 
of the Bridge. The Model 2 scenario predicts that the ma- 
jority of the Bridge material originated from the SMC, but 
there should be some contribution from LMC gas that in- 
creases in importance with proximity to the LMC. Model 2 
thus predicts that Bridge material towards the LMC should 
be increasingly metal enriched. This should not be true in 
Model 1. 

In the edge-on view of the Model 1 disk, the gas disk ap- 
pears to be tilted relative to the stellar distribution. This is 
likely because of the infall of gas from the SMC that forms 
the Magellanic Bridge. The outer stellar and gas disks in 
Model 2 are significantly warped and distorted in the edge- 
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Figure 10. The line-of-sight distances for the gas distribution of the simulated system are plotted as a function of Magellanic Longitude 
for Model 1 (left) and Model 2 (right). The modeled line-of-sight distance along the orbit of t he L MC(SMC) is plotted as the solid 
(dashed) yellow line. The solid red line indicates the distance estimate from jjin Lvnden-Bell (|200^ \ 



Figure 11. The stellar surface density of the LMC disk is plotted in the line-of-sight frame for Model 1 (top) and Model 2 (bottom). 
RA and DEC grids are overplotted in green across the face of the disk. The stellar distribution in Model 2 is significantly more disturbed 
than in Model 1. In particular, the bar in Model 2 is off-center relative to the stellar disk. The Model 2 LMC disk is also significantly 
disturbed in the upper left (North-East). 



on view. Th e true disk is also obse rved to be both flared 
(lAlves &: Nels on 20 0d) and warped (Ivan der Marel Cionil 
I2OOII : lQlsen"fc Salvy[200i : iNikolaev et al.""2004). Such re- 
su lts are in keeping wit h a study of Magellanic-type spirals 
bv lWilcots et al.1 (| 19961 ). who also suggest that the observed 
lopsidedness in their HI disks may be a result of minor merg- 
ers. 

In Model 2 the gas disk has also formed a pronounced 



arc in the upper right. Since our star formation prescrip- 
tions depend sensitively on the gas density, this arc of gas 
will also be actively forming stars (see ^ 16. ip . giving the LMC 
the appearance of a one-armed spiral. A number of numer- 
ical studies have been conducted on the resulting structure 
of a large galaxy after a direct collision with a smaller com- 
pani on in the context of explaining the origin of rin g galax- 
ies ((Lvnds Toomrelll976l : IWeil &: Hernauistlll993l : IStruckl 
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Figure 12. The gas (left) and stellar (right) density of the LMC disk for Model 1. Top: face-on view (x,y plane). Middle: edge-on view 
along the x axis. Bottom: along the y axis. Both the gas and stellar projections have a centered, in-plane bar. 



Il997l l. In particular, IStruckl (Il997l ) finds that, in some cases, 
a one-armed spiral structure can be excited in the larger 
galaxy. Also, Bekki (200 9) explored a scenario where the 
LMC bar becomes off-center as a result of a recent encounter 
with a dark 10^ companion. The specific asymmetries in- 
duced depend sensitively on the mass ratio, inclination and 
the location of the smaller companion's passage through the 
larger galaxy. A future study will explore these parameters 
in depth in the context of the LMC-SMC encounter and 
assess the lon gevity o f the resulting asymmetric structures. 
For example, iLevine Sparke (J998,) have illustrated that 



disk-lopsidedness can be long-lived if the disk is displaced 
from the center of the dark matter potential and spinning 
in a sense that is retrograde to its orbit about that center. 



Finally we note that, while we have not discussed the re- 
sults for the SMC in detail, the simulated SMC morphology 
in Model 2 is consistent with the observations of a "bar" -like 
main body with a stellar wing leading towards the LMC and 
a significant line-of-sight depth. We will discuss these results 
in depth in a future paper (Besla et al. 2012 in prep.). 
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Figure 13. Same as Figure [12] but for Model 2. The stellar bar is clearly off-centered and warped relative to the stellar disk; it is also 
absent in the gas. The SMC's orbit in an LMC-centric frame is overplotted (dashed line) in the middle and bottom panels of the LMC's 
gas distribution. Various times are also marked along the SMC's past orbit. The SMC collides with the LMC 100 Myr ago and the SMC's 
current position is marked. An extra-planar stream of gas is pulled out from the LMC by the passage of the SMC. 



5 LMC KINEMATICS 

The internal kinematics of the LMC has been quantified 
by many tracers. The LMC's rotation curve rises roughly 
linearly to a radius of ^ 4 kpc, after which it stays flat 
at a value of Vrot- The observed rotation curve has been 
noted to peak at different values depending on the kine- 
matic tracer bein^ studied. The HI kinematics yield Vrot = 
80 km/s ( St avelev- Smith et al.l [2003a ), data from red su- 



pergiants gives Vrot = 107 km/s (lOlsen &: Massev 2007 ) 
and carbon stars yields V rot=61 km/s (Ivan der Marel et al.l 
l2002h . However, recently lOlsen et al.l (|201lh examined the 
kinematics of a combined population of massive red super- 
giants, oxygen-rich and carbon-rich AGB stars in the LMC. 
After correcting for the LMC's space motion and the asym- 
metric drift in the AGB population, they find a consistent 
rotation curve between all kinematic tracers with Vrot = 87 
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±5 km/s. This is in accord with the HI rotation curve and 
the initial value adopted in this study 95 km/s). 

The LMC disk initial conditions (gas fraction, equation 
of state) are chosen such that the LMC kinematics are rep- 
resentative of a symmetric, bar-unstable disk galaxy. It re- 
mains to be seen whether MW tides will introduce kinematic 
anomalies in the disk in a first infall scenario, or, perhaps 
more significantly, whether the LMC can retain a kinemat- 
ically stable disk after a direct collision with the SMC (i.e. 
in Model 2). 

In Figures 1141 and [15] the kinematics of the LMC disk 
in Models 1 and 2, respectively, are broken down for various 
kinematic tracers: gas (top panel), young stars (middle) and 
old stars (bottom), in the line-of-sight frame. In each panel 
the disk is centered on the stellar center of mass. The cen- 
ter panels show the surface density of the tracer population. 
A slit is placed along the largest velocity gradient of the 
LMC's older stellar distribution and is defined as the major 
kinematic axis (red). A second slit is placed 90 degrees with 
respect to the major axis and referred to as the "minor" 
kinematic axis (blue). The position angle of the kinematic 
major axis of the simulated disk for Model 1 is 55 degrees 
counter-clockwise from the x-axis in all panels and 50 de- 
grees for Model 2. The line-of-sight velocities along the slits 
are plotted in the left panel. The right panel shows the full 
line-of-sight velocity field. The middle and right boxes are 
18 kpc a side whereas the left box is scaled to the length of 
the slit (16 kpc along the x-axis). 

In both models the LMC disk retains a pronounced ve- 
locity gradient along the same major axis; perhaps surpris- 
ingly, the LMC disk retains a well defined rot ation curve de- 
spite a direc t collision with the SMC. Indeed, [Hopkins et al.l 
(|2008l . 12009') showed that disks can survive even a 1:1 mass 
ratio major merger. However, the disk kinematics in Model 
2 are more distorted than in Model 1, particularly the 
velocity field. 

There are observed asymmetries in the LMC's gas and 
stellar kinematics. It has long been noted that the HI kine- 
matic center is offset by ^1 kpc from both the stellar kine- 
matic and photometric center, which is roughly c entered on 
the stellar bar (as illustrated in ICole et al.ll2005l ). However, 
upcoming work by Kallivayalil et al. (in prep) using a 3rd 
epoch of HST data provides proper motions of high enough 
accuracy to independently constrain all parameters of the 
LMC rotation field and geometry, including the dynamical 
center. The best-fit stellar dynamical center from the proper 
motions agrees with the HI dynamical center determined by 
Kim et al. (1998), but remains offset from the photometric 
center (van der Marel & Kallivayalil, in prep.). 

Each panel in Figures [TJ] and [15] is centered on the 
peak of the stellar density of the simulated LMC disk; i.e. the 
photometric center. In Model 1, the stellar density peak is 
coincident with the kinematic centers of all tracers (vertical 
dashed line in left panel crosses zero where the major slit 
velocities do). Thus, contrary to observations, the stellar and 
gas kinematic centers are coincident with the photometric 
center and the center of the stellar bar. 

In Model 2, the zero velocity field of the stars and gas is 
twisted such that the velocity gradient does not cross zero 
at the location of the stellar density peak. The kinematic 
centers of all tracers are offset by about 1 kpc from the pho- 
tometric centers, as observed. However, the gas and stellar 



kinematics are somewhat discrepant from each other. This 
is illustrated by placing a third slit 1 kpc above the ma- 
jor axis (illustrated in the right hand panel of Figure 1151 
in green). The old and young stellar line-of-sight velocities 
along this slit are similar to those along the original major 
axis. The gas line-of-sight velocities, on the other hand, have 
changed, crossing the zero axis roughly 1 kpc further away 
than seen for the stars. This implies that the shape of the 
zero velocity field across the face of the gaseous and stellar 
disks are different and is likely a result of the warped stellar 
bar. Note that neither model predicts strong differences in 
the rotation curves traced by the you ng (< 1 Cyr) and o lder 
stellar populations, as expected from Olsen et al.l (|201lh . 

Again, the simulated SMC kinematics will be presented 
in a forthcoming paper. However, we mention here that the 
resulting kinematics are much more consistent with Model 
2 than Model 1. It appears that a direct collision with the 
LMC is required to erase the initial velo c ity gr adient in the 
older stellar population. IZaritskv et al ] (|2000h was unable 
to find a pronounced velocity gradient in the RGB popu- 
lation within a radius of 2 kpc of the center of mass, de- 
spite the existence of a 6 km/s velocity gradient in the gas 
(|Stanimirovic et al1l2004 V In our simulations, the gas is able 
to cool, since it is dissipative, and the original gas disk sur- 
vives the tidal shocks resulting from the direct impact with 
the LMC, whereas the older stellar population does not. In 
Model 1 (no direct impact with the LMC), there is always 
a pronounced gradient in both the stellar and gaseous com- 
ponents of the SMC. 



6 DISCUSSION 

In this study we have shown that it is possible to explain 
the nature of the LMC's off-center stellar bar (§|4|) and gas 
and stellar kinematics (§[5]) in a model that self-consistently 
reproduces the general large scale gas morphology of the 
Magellanic System (§[3]) in a first infall scenario. To do this, 
we invoked a recent direct collision between the LMC and 
SMC (Model 2). Here we discuss some of the testable con- 
sequences of such a scenario. 

Note that the resulting simulated SMC kinematics and 
structure and an expected stellar counterpart to the Magel- 
lanic Stream will be discussed in upcoming papers. 

6.1 The Recent Star Formation History of MCs 

A direct recent collision between the LMC and SMC 
would likely leave notable marks in the star formation 
histories (SFHs) of both of these galaxies. A correlated 
burst of star formation during such a recent encounter has 
been theorized in many previous numerical studies (e.g., 
IGardiner Noguchilll996l : iBekki Chiball2005l . l2007h . Ob- 
servationally, there is significant debate over the existence 
of correlated bursts of st ar formation within both galaxies. 
iHarris Z aritskv (2009) claim that the total SFR in the 
LMC was higher than average ^100 and 500 Myr ago and in 
the SMC at ^60 and 400 Myr ago. Other authors claim that 
the LMC shows global enhancements ~ 125 and 800 Myr 
ago (jPietrzvnski Udal ski 2000al ) and i n the SMC at -100 
Myr only fPietrzvn ski Ud alski 2000b'). The general con- 
sensus does appear to be that both galaxies show a steadily 
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Figure 14. The kinematics of the simulated LMC's gaseous disk (top), young stellar disk (< 1 Gyr; middle) and old stellar disk (> 
1 Gyr; bottom) in Model 1 are illustrated in the line-of-sight frame (North to the top, East to the left; the SMC is located to the 
South- West). The right panels show the line-of-sight velocity field, where the center-of-mass velocity of the respective kinematic tracer 
has been subtracted. Density contours and major axis slit location are also indicated. The color gradient denotes material moving towards 
(blue) and away (red) from the observer. The central panels show the surface density of the tracer, with density contours overplotted. 
The central and right panel boxes spans 18 kpc a side. Each box is centered on the stellar density peak. Slits are placed along the major 
(red) and minor (blue) stellar kinematic axes as indicated. The major kinematic axis slit is inclined 55 degrees counter clock wise from 
the X-axis in all panels. The minor kinematic axis is placed 90 degrees with respect to the major axis, although in practice this does not 
exactly trace the zero velocity regions. Note that even if the gas kinematics are being examined, the red slits still denote the locations of 
the stellar kinematic axes. In the left panel, the line-of-sight velocities are plotted along the slit. The dashed vertical line indicates the 
location of the peak stellar density (photometric center); in Model 1 it is coincident with the kinematic centers of all tracers. 



increasing; SFR o yer the past Gyr 



iNoel et alJlixil : iMcCumber et al.l 



Harris ZaritskYll2009l : 



2005l l. Testing this sce- 



nario reliably depends on the accuracy of the timing of the 
collision/interactions between the L/SMC and the adopted 
star formation prescriptions. 

We can bracket the time range for this impact to 
be within 100-300 Myr. The upper limit on the collision 
timescale comes from the oldest detected stellar populations 
i n the M agellanic Bridge, which are believed to form in-situ 
(|Harrisl 2007) and thus mark the formation of time of the 



Bridge. As such, given the model parameter uncertainties, 
we cannot use the models here to definitively predict the 
timing of the collision and consequent star formation. 

The exact timing of the collision is strongly dependent 
on model parameters. In the presented Model 2, the most 
recent collision occurred ^100 Myr ago. However, the MCs 
are too close together today in this model, indicating that 
the true collision likely occurred earlier. 

Given our adopted prescriptions, the modeled SFHs of 
both MGs are plotted as a function of time since they first 
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Figure 15. The same as in Figure [M] except for Model 2. The kinematic major axis (red) is inclined 50 degrees counter clock wise 
from the x-axis in all panels. In the left panel it is clear that the kinematic center of all tracer populations are coincident, but not with 
the photometric center (dashed vertical line). Furthermore, the shape of the velocity field as traced by gas vs. stars is different across 
the face of the disk. This is illustrated by the inclusion of a new slit (in green), placed 1 kpc above the major axis. The line-of-sight 
velocities of the gas along the 'New' slit deviate from their stellar counterpart. This offset is likely a result of the warped stellar bar. 



Table 3. Current Star Formation Rate 



Galaxy 


Model 1 


Model 2 


Ra & IR 


Free-free Emission 




(Mo/yr) 


(Mo/yr) 


(Mo/yr) 


(Mo/yr) 


LMC 


0.5 


0.4 


0.25 


0.14 


SMC 


0.2 


0.4 


0.05-0.08 


0.015 



Note. — The SFR is computed within a radius of 15 kpc for the LMC and 2 kpc for the SMC. The third co lumn indica t es the 
mean values determined from Ho: and MIPS emission bv lWhitnev et al.l (|2008i ) for the LMC and from IR and Ha bv lWilke et al.l (|2004j ) 
and Kennicutt & Hod gj (Il986l ), respectively. The fourth column presents a lower limit, determined from the free- free flux measured by 
iMurrav RahmanI (j201ol ) using WMAP. 
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Figure 16. The simulated SFR is computed within a radius of 15 kpc for the LMC (sohd hues) and 2 kpc for the SMC (dashed hues) and 
plotted as a function of time since the MCs first crossed within R200 of the MW for Model 1 (red) and Model 2 (blue). The corresponding 
orbital histories are plotted in Figure |4] In Model 2, the MCs experience a close encounter at 0.1 Gyr ago (time indicates today). 



crossed within R200 of the MW in Figure [161 Note that the 
star formation rate is derived from the gas density above a 
set threshold value. Contrary to Model 1, the SFH of the 
SMC o ver the past Gyr in Model 2 increases ste adily, as ob- 
served ("Noel et al. 2QQ9l: lMcCumber et al.lliooBI ). The Model 
2 result occurs because the separation between the MCs is 
smaller than in Model 1, hence the relative importance of 
tidal distortions from the LMC and consequent triggered 
star formation is also stronger. 

The LMC is observed to be unusually blue relative 
to ana logs with similar R -band magnitudes identified in 
SDSS (i Tollerud etHI l201lh . This is likely a result of two 
factors: 1) if the MCs are on their first infall and just 
past their first pericentric approach to our MW, they may 
be experiencing triggered star formation induced by MW 
tides. 2) Interactions between the MCs have likely kept 
the SFR in the LMC higher than it would be if it did 
not have a companion. The SFRs of galaxies are known 
to i ncrease as a function of separation to a close compan- 
ion (iLarson TinslevI Il978l : iFreedman Woods et all |201Q| : 
IPatton et al.ll201lL and references therein). Since L MC-SMC 
pairs are rarely found around MW type hosts (|Liu et al.l 
I2OIII : iBovlan-Kolchin et al.l 120111 ). it is natural that the 
LMC should have an anomalously high current SFR rela- 
tive to the average analog (which is more likely to be iso- 
lated). This theory can be tested by comparing the LMC's 
color/SFR to a sample of Magellanic Irregulars (LMC 
analogs) with known close companions. 

The recent collision between the MCs in Model 2 re- 
sults in a dramatic increase in the SFR in the SMC at that 
time. As discussed above, the exact timing of the true col- 
lision is quite uncertain. The magnitude of the simulated 
burst is inconsistent with obs ervations of the SFH of the 
SMC (jHarris k Zaritskvl l2006i ) . Furthermore, the modeled 
SFRs today are also higher in both MCs than observed (see 
Table [3|) . This likely points to significant problems in our 
adopted star formation prescriptions. 

The star formation prescription we have adopted in 
these simulations follows a Kennicutt-Schmidt volume den- 
sity law with a local volume density cutoff for star for- 



mation of uh — 0.13cm' ISpringel k HernauistI (I2OO3I) 
showed that when combined with appropriate star forma- 
tion timescales and typical scale-heights this gives a good 
match to the observed Kennicutt star formation surface den- 
sity relation for relatively massive galaxies. However in the 
regime where galactic gas is dominated by atomic hydro- 
gen and where molecular, star forming, gas constitutes only 
a fraction of the gaseous content, one has to properly ac- 
count for the formation of loc al density enhancements and 
molecular hydrogen f ormation (iRobertson k Kravtsovll2008l: 
Kuhlen et all I2OIII : iKrumholz et al.l l2009l : iHopkins et al.l 
20 111 ). In low mass galaxies th at tend to be metal poor, such 
as the SMC (|Fox et al.l l2010'), molecular gas formation and 
therefore also the star formation is very inefficient. Theoreti- 
cal models show that not accounting for the details of molec- 
ular gas formation and using global metallicity independent, 
low density, threshold for the ISM model can lead to seri- 
ous overestima tes of th e star formation rates of metal-poor 
galaxies (Kuh len et al Eolll While larger improvement can 
be made by accounting for metallicity dependent molecular 
gas formation in sub-resolution models, direct modeling of 
processes that self-regulate formation of molecular clouds, 
their star formation and related feedback requires more a 
complex ISM model and re solution that is beyo nd the sim- 
ulations used in this work ([Hopkins et al.ll201ll ). 

Since the bulk of this study focuses on the morphologi- 
cal and kinematic properties of the simulated MCs and the 
dynamics of the Magellanic System, the detailed star forma- 
tion prescription does not alter any of the main conclusions 
in this work. However, it does limit the predictive power of 
the models concerning the chemical evolution and SFHs of 
the MCs. Future detailed studies of star formation prescrip- 
tions/feedback in a repeated series of Model 2 collisions can 
be compared directly with the multiples observational data 
sets for the SFHs of the MCs. Such a study may be a pow- 
erful method of constraining the appropriate sub-resolution 
physics for Magellanic Irregulars and shed light on to the 
nature of bursts of star formation in the histories of the 
MCs. 
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6.2 Star Formation in the Bridge and Stream 

An additional check for our proposed Models is the pre- 
dicted locations of ongoing star formation in the simulated 
Magellanic Bridge and Stream. 

Figure [TTl shows the instantaneous SFR density in the 
Magellanic System for Models 1 and 2; the SFR density is 
derived from the gas density. These plots indicate the loca- 
tion of gas with densities above the star formation threshold 
and should thus be forming stars. While these plots may not 
be quantitatively accurate (the results clearly depend on the 
star formation prescription), they do highlight the location 
of high density gas that should be forming stars in any sub- 
grid model. 

In both cases the gas densities in the Magellanic Stream 
are too low to form stars. This is consistent with observa- 
tions; molecular gas has not been detected in the densest 
cloudlets of the Stream, sug gesting that star form ation is 
not actively occurring there ([Matthews et al.ll2009l V 

In Model 1 stars do not form in the Magellanic Bridge, 
whereas in Model 2 a well-defined bridge of star forming gas 
is seen connecting both galaxies. Star s as young as 10-40 
My r have been detected in the Bridge (.Demers Battinellil 
1 19981 ), indicating that star formation is on-going there, as 
these stars could not have migrated from the SMC dur- 
ing their lifetimes. Furthermore, Harris (2007) was unable 
to locate stars older than 300 Myr in the Bridge, lending 
support to the idea that the majority of young stars there 
were formed in-situ. The different results between Model 1 
and Model 2 indicate that a very close/direct encounter be- 
tween the MCs may be required in order to trigger in-situ 
star formation in the Bridge, otherwise gas densities in the 
tidal bridge should not be significantly different from the 
Magellanic Stream (i.e. the tidal tail). A close encounter 
has generally been invoked to explain the prope rties of the 
Bridge in previous st udies (e.g., lYoshizawa No guchi 2003; 
IConnors et al.ll200^ . Clearly shock induced star formation 
is a relevant process in the direct collision scenario presented 
here and must be accounted for in order to characterize this 
environment accurately. 

Note that, although these results seem to favor Model 2, 
ram pressure compression of the Bridge region may also lea d 
to triggered star formation fe.g.. Mastropietro e t al.ll2009 V 
Although, in this case there likely should be stars forming in 
both the Leading Arm and the Magellanic Stream as well. 

The SFR along the Bridge in Model 2 increases steadily 
towards the SMC; th is should be testable observationally. 
Indeed, iHarri 

i 12007) showed that the distribution of blue 
stars is more dense in the Bridge towards the SMC. Spitzer 
observations also confirm the presence of you ng stars in the 
"wing" of the SMC, which leads to the Bridge ([Gordon et al.l 
I2OO9I V Furthermore, Ra measurements of the Bridge by the 
Wisconsin Ha Mapper (WHAM) survey indicate that the 
Ha emission in the Bridge increases steadily towards the 
SMC (Barger et al. in prep). 

In Figure [17] the bar region of Model 1 is currently the 
dominant site of ongoing star formation in the LMC (con- 
trary to observations) . The distribution of star formation in 
the disk appears similar to an isolated galaxy with a num- 
ber of spiral arms. On the other hand, in Model 2 the bar 
is not the most prominent site of star formation today. Fig- 
ure [18] provides a zoomed in view of the Model 2 results 



in the line-of-sight frame for the LMC disk. The center of 
mass of the stellar disk is located just to the right of an in- 
tense star-forming knot that indicates the impact location 
of the SMC-LMC collision. As discussed earlier, this occurs 
because the bar is warped with respect to the gas disk plane 
and is therefore inefficient at funneling gas along it. 

Without appropriate feedback prescriptions the appear- 
ance of the impact site today is difficult to assess. The thick- 
ness of the gas disk near the impact site has increased rela- 
tive to the rest of the disk (see where the SMC's trajectory 
crosses the LMC disk plane in Figure [13]) . Intriguingly, this 
i s also true of the active 30 Doradus star forming region 
([Padoan et al.llioOll V which is located at roughly the same 
location. However, the SFH of the true 30 Doradus region 
indicates that it has only been an active site of star forma- 
tion for the past 12 Myr ([Harris k Zaritskv 2009). Without 
accurately capturing star formation induced by shocks and 
including feedback, it is unclear whether the initial collision 
(100-300 Myr ago) is related to the triggering of star forma- 
tion in 30 Doradus today or whether its remnan t is related 
to an y of the sup ergiant shells in the LMC disk ([Kim et al.l 
fl998l : iBook et al.ll2008l . l2009l ) . 

6.3 Mass Breakdown 

The initial and final mass contained within characteristic 
radii of the simulated LMC and SMC and the resulting mass 
of the simulated stream are listed in Table [J] and compared 
with observations. Quoted gas mass estimates are for the 
total gas component; the neutral HI content will be lower 
than these values. We do not include a UV ionization model 
in these simulations and thus cannot accurately estimate the 
neutral gas fraction. 

In general, the mass estimates for the LMC and SMC 
agree with the observations within a factor of two. However, 
the final gas mass estimates for the MCs are lower than those 
observed and the simulated stream is about a factor of five 
lower. 

As discussed in § 16.11 given the adopted star formation 
prescriptions, the SFR in the MCs is being overestimated in 
the simulations. As a result, the gas consumption timescale 
is also overestimated. The adoption of different prescriptions 
and/or changes in model parameters so that the SMC does 
not lose as much gas at early times (e.g. by reducing the 
scale length of the SMC's gas disk) may allow for the SMC 
to retain its gas for longer, allowing more material to be 
removed at later times. 

The mass budget of the Stream is also likely underesti- 
mated because ram pressure stripping has not been modeled 
here; the bulk of the missing material is in the region clos- 
est to the MCs, rather than the tip of the tail. But this still 
means that more gas needs to be retained by the MCs them- 
selves to provide a gas reservoir for this process to operate 
at late times. This suggests that the larger problem here is 
the adopted star formation prescriptions. 

Probably more significantly, the total initial gas budget 
of the MCs has been underestimated in these simulations; 
the amount of gas initially modeled in both the LMC and 
SMC cannot explain the total gas budget (including both 
neutral and ionized components) of the Magellanic System. 
Although the neutral HI gas content of the observed Magel- 
lanic System has been well quantified, the ionized gas frac- 
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Figure 17. The instantaneous SFR density for the simulated Magellanic System is plotted for Model 1 (top panel) and Model 2 (bottom) 
as a function of Magellanic Longitude. The SFR density is proportional to UV luminosity and identifies the location of gas that is currently 
forming stars. The white circle indicates the observed extent of the LMC's disk and the solid white line denotes the observed location of 
the Stream. Blue solid(dashed) lines indicate the past orbits of center of mass of the LMC(SMC). 




Figure 18. A zoom-in on the instantaneous SFR density in the LMC and SMC in Model 2 from the bottom panel of Figure fTTl RA 
and DEC grids are overplotted in green (North is to the top and East is to the left). This image indicates where stars are expected to 
be forming in the LMC today: star forming gas is concentrated in the South-East and a single spiral arm of star formation is seen in the 
North- West. The SFR increases along the simulated bridge towards the SMC (South-East). 
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tion is poorly understood. IFox et al find that at the 

tip of the Stream as n iuch as 95% o f the gas may be ionized. 
Along a similar line, i Lehner et al find that gas in 

the Magellanic Bridge is 70-90 % ionized. This means that 
the total gas mass budget of the Stream (and consequently 
of its progenitor) is significantly underestimated. Note that 
this statement implies that the LMC and SMC must have 
a significant amount of dark matter in order to make their 
initial baryon fractions consistent with cosmological expec- 
tations; this is the main motivation for the large total infall 
masses adopted in this study. 

The fact that the stellar mass in the central regions 
of the SMC does not change substantially indicates that 
the bulk of the material being removed from the SMC via 
LMC tides is from the outskirts of the SMC. This explains 
why the stellar counterpart of the Stream is so faint. We 
will comment on the observability of the this stellar stream 
counterpart in an upcoming paper. 

6.4 Tidally Captured SMC Stars 

In the presented models there is a continual transfer of ma- 
terial from the SMC to the LMC. In particular, there is ex- 
pected to be a population of stars that are tidally stripped 
from the SMC and captured by t he LMC in both models. 

Recently, lOlsen et al.l (|201lh discovered a population 
of metal poor RGB stars in the LMC field that have dif- 
ferent kinematics from those of local stars in the LMC 
disk. [Craff et al. (2000) have also identified a possible kine- 
matically distinct collection of carbon stars and suggest 
that this population lies outside of the LMC. The dis- 
covery of such stars is a natural theoretical expectation 
from any tidal model for the LMC-SMC interaction. To 
date, no stars have been detect ed in the Magellanic Stream 
(jCuhathakurta Reit^ll998l l and only stars younger tha n 
300 Myr have been id e ntified in the Bridge (|Harrisl lioOTh . 
Note that the iHarri i (|2007l l observations focused on the 
leading ridgeline (location of the highest gas density) of the 
Magellanic Bridge. This leading edge would currently be ex- 
periencing maximal ram pressure and so it is possible that 
the peak gas density in the bridge is displ aced fr o m the tidal 
stellar population theorized to be there. iHarrisI (l2007l l con- 
strained the stellar density of a possible offset stellar popula- 
tion using 2MASS observations, but the 2MASS sensitivity 
limit of 20 Ks mag/arcsec^ is likely far too low to detect the 
expected faint stellar bridge (>30 Ks mag/arcsec^; Besla et 
al. in prep). 

The potential identification of tidally stripped SMC 
stars in orbit about the LMC may be a key discriminant 
between various model interpretations of the origin of the 
Magellanic Stream, as they should not be present in a pure 
hydrodynamic model (Mastropietro et al. 2005) or one that 
relies on stellar outflows (Nid ever et al.l l2008l : lOland 12004 ). 

To test whether the sim ulation resul t s for the tidal de- 
bris are consistent with the lOlsen et al.l (|201lh detections, 
we plot in Figure [19] the expected distribution and kinemat- 
ics of the stars captured by the LMC from the SMC for both 
Model 1 (top panel) and Model 2 (bottom panel). The box 
size and orientation is the same as in Figures [14] and 1151 
that is, the field of view is centered on the LMC. The stel- 
lar line-of-sight velocities have been corrected for the center 
of mass motion of the LMC. The SMC is located towards 



the South- West in this viewing perspective. Note that the 
SMC is actually present in this field-of-view for the Model 2 
results, as this simulation resulted in the MCs being closer 
together than observed. 

In Model 1 these transferred stars form a well-defined 
arc that is in orbit about the LMC: these stars are located 
behind the LMC disk. Comparing the velocity field to the 
bottom-right panel of Figure I14[ the kinematics of these 
stars appear to be offset by nearly 90 degrees relative to 
the velocity gradient in the LMC's disk. In Model 2, the 
stellar debris from the SMC exhibits a large range in veloc- 
ities (± 150 km/s). We made a velocity floor and ceiling of 
±70 km/s to better compare to the LMC stellar disk kine- 
matics. In the North- West there are stars that appear to be 
moving towards the observer and in the South-East there 
are stars moving away from the observer - this is opposite 
to the observed kinematics of the LMC stellar disk. Thus in 
both models, stellar debris captured by the LMC from the 
SMC is expected to have kinematics that are distinct from 
those of the LMC disk stars - this is a generic prediction of 
the BIO model and is consistent with the Olsen et al. (201l|) 
obs ervations. 

lOlsen et al ] (|201lh estimate the mass of the observed 
SMC debris population to be 5% of that of the LMC's cur- 
rent disk mass (i.e. ^ 1.4 x 10^ M©). In Model 1, only a 
modest amount of stars is transferred from the SMC to the 
LMC (^ 0.2%). As such. Model 1 cannot account for this 
population. On the other hand, in Model 2, the LMC ac- 
cretes 1.5% of its current disk mass from the SMC. 

The exact distribution of SMC debris is certainly depen- 
dent on a number of parani eters, but overall we can conclude 
that the lOlsen et al ] (|201lh results support a model in which 
LMC tides have been actively distorting the SMC. Further- 
more, in both models the accreted stars are largely located 
behind the LMC disk and may provide a natural explanation 
for the origin of the observed MACHO microlensing events 
(Besla et al. in prep). 

6.5 The Nature of Magellanic Irregulars 

Dwarf galaxies are broadly referred to as galaxies with lumi- 
nosities <0. 1-0.3 of L*. This definition encompasses a wide 
range of objects of varying morphology, including both the 
LMC and SMC. In this work we have introduced a mass 
model for the LMC with a total mass of ^ 10^^ M©; it is 
questionable as to whether such a massive galaxy should 
be included in the same category as dwarf Spheroidal or 
dwarf Irregular s galaxies. In particular, deta iled analysis of 
the geometrv (Ivan der Marel Cionil 1200 ih and kinemat- 
ics (jvan der Marel et al.ll2002l ) of the LMC prove that it is 
a disk galaxy. Furthermore, when looking at the distribu- 
tion of intermediate-age and old stars out to large radii (i.e. 
ignoring the visible light in the bar region), the LMC does 
not look at all irregular, but clearly res embles a spiral galax y 
wit h an asymmetric one-armed s piral (|van der Marelll2001h . 

'de Vaucouleurs Freeman] (|l972l ) suggested that Mag- 
ellanic Irregulars have more in common with spiral galaxies 
than dwarfs, referring to the LMC as an asymmetric, late- 
type barred spiral galaxy. They argue that Magellanic Irreg- 
ulars represent an extension of the Hubble spiral sequence 
(Sa, Sb, Sc, Sd, Sm, Im), where the subscript m denotes 
"Magellanic" (see also, section 4.1.1 of iBinnev &: Me rr ifieldl 
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Table 4. Initial and Final Mass Distributions 



vjj cxi cxjv y 


P T O p y i" A 7" 


In-itjiOjl 


Model 1 


Model 2 


L^oC-i V C-vJ. 






(-1 n9 ]wr _ ^ 


yiu ivi0; 


ivi0j 




LMC 


Stars (< 9 kpc) 


2.5 


3.1 


3.1 


2.7 (1) 




Gas (< 5 kpc) 


0.87 


0.17 


0.26 


0.441 (2) 




Total (<9 kpc) 


18 


21 


21 


13(±3) (3) 




Total (<4 kpc) 


6.9 


7.3 


7.1 


5 (4) 


SMC 


Stars (< 3 kpc) 


0.20 


0.26 


0.20 


0.31 (5) 




Gas (< 3 kpc) 


0.20 


0.16 


0.18 


0.42 (5) 




Total (< 3 kpc) 


2.1 


1.8 


1.1 


2.7-5.1 (6) 












2.4 (4) 




Total(< 1.6 kpc) 


0.84 


0.70 


0.49 


1.4- 1.9 (6) 


Stream^ 


Gas 




0.12 


0.10 


0.5 (7) 


Stars 




0.6 


0.4 





^ The Stream is defined as material at Magellanic Longitude less than -30. We have also accounted for the distance in our 
simulated stream in order to properly compare to the observed mass, which is computed modulo the distance squared. 



Note. — (1) Ivan der Marel et aP (l2002l) : The out e rmost data point lies at 9 kpc;_(2) iBriins et aP (|2005l ): (3) Ivan der Marel et al.l 
(|20Q9l ) ; f4) lKim et al.l (1998): r~ IStanimirovic et al.l hOOd ): (6) [karris & Zaritskv (20 oH); (7) Her e we have accounted f or the average 
dista nce to the simulated streams. The measured value is 1.25 x 10^ ((i/55kpc)2 (Br iins et al.ll2005l ) + 2.0 x 10'^(d/120)2 (|Nidever et al.l 
I2OIOI V The quoted simulated gas mass refers to the total mass, not just that of the neutral HI, whereas the observed values refer to HI 
gas only. 




R (kpc) R (kpc) 

Figure 19. The kinematics of SMC stars that are captured by the LMC in Model 1 (top panel) and Model 2 (bottom panel) are 
illustrated in the line-of-sight frame of the LMC's disk. Only particles initially belonging to the SMC are plotted. The center-of-mass 
velocity of the LMC's stars has been subtracted from the velocity field of the transferred material. The left panels shows the stellar 
density map with the highest density contours overplotted. The size of each box is 11 kpc per side, the same as that in Figures 1141 and 
1151 . The slit along the major and minor kinematic axes of the LMC are overplotted for reference. The color gradient ranging from blue 
to red indicates material moving towards (blue) and away (red) from the observer. 
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1 19981 ). Spiral structure "decays" along the sequence, with Sc 
having irregular spiral patterns and Im none at all. This is 
true also of the Barred Spiral sequence (SB), with asym- 
metry referring also to the appearance/location of the bar. 
Magellanic Irregulars encompass the late stages of both 
barred and unbarred spirals (e.g. Sd-Im, SBd-SBm); the 
LMC is classified as SBm under this scheme and the SMC 
an Im. 

This work postulates that many barred Magellanic Ir- 
regulars may be perturbed versions of symmetric low-mass, 
bulgeless, barred galaxies, such as SBc type galaxies, where 
the bar is typically well- centered. We further illustrate a 
mode of inducing such perturbations, namely interactions 
with lower mass companions. In this picture, the LMC 
should not be thought of as a "dwarf galaxy except in the 
sense that it is less luminous than the MW. 

Asymmetric bars are typically not seen in massive 
galaxies, making them a defining characteristic of Magel- 
lanic Irregulars. This might be explained by noting that 
Sc/SBc type galaxies and Magellanic Irregulars (Sd-Im) are 
low mass systems compared to MW type galaxies and, cor- 
respondingly, they sample very different environments. Such 
galaxies do not have bulges, have shallower central poten- 
tials, are more dark matter dominated, have different halo 
concentration parameters, higher gas fractions and different 
ratios between the ISM temperature and virial temperature 
of their halos: all of these differences will influence the re- 
sponse of the system to tidal perturbations. As such, even 
if the physical scenario is similar (mass ratio and orbital 
configuration), the response of a MW type galaxy to a 1:10 
mass ratio direct collision is expected to be different than 
an LMC-SMC (^1:10 mass ratio) encounter. 

It is possible, for example, that the presence of a bulge 
may aid in stabilizing the bar of a high mass galaxy, prevent- 
ing comparable asymmetries in the bar from arising. How- 
ever, a detailed numerical study of such parameter space to 
test such conjectures is beyond the scope of the presented 
study. 

Magellanic Irregulars are ubiquitous in our Local Vol- 
ume. In light of the theory presented here, it must also 
be true that interactions between low-mass barred galax- 
ies and sma ller companions are a relatively frequent occur- 
rence. From Hopkins et al. ](|2010^, the expected galaxy ma- 
jor merger (mass ratio> 1/3) rate is relatively flat as a func- 
tion of host galaxy mass at z=0, so it is not expected that 
such encounte rs would be more likely for low mass systems. 
IStewart et all (|2008) find that 25-40% of hosts with mass 
of order ^ 10^^ Mo/h have accreted a 1:10 mass ratio sub- 
halo within the past 6-8 Gyr, which is the timescale for our 
isolated LMC-SMC interaction. This scenario is thus not at 
odds with cosmological expectat ions. 

It has been pointed out by IWilcots PrescottI (l2004l l 
that many Magellanic Irregulars do not currently have com- 
pani ons, despit e earli er claims of a high frequency of pairs 
by lOdewahnI (|l994l ). However, the number of observed 
interacting dwarf systems is s teadily increasing. Recently, 
iMartinez-Delgado et al.l ( j201 1) have discovered a stellar 
stream about the Magellanic Irregular galaxy NGC 4449, 
which is an LMC analog in terms of its absolute magnitude. 
Although the stellar mass ratio of the disrupted object and 
the host is 1:50, the inferred dynamical mass ratio is between 
1:10 - 1:5, making this system an analog of the late stages 



of an LMC-SMC type tidal interaction. NGC 4449 was long 
thought to be an isolated Magellanic Irregular until observa- 
tions of associated HI streams indicated that it likely had a n 
encounter with an unseen companion (jHunter et al1ll998l ). 
Moreover, an unusual globular cluster exists in this galaxy 
with pr operties consistent wi th the nucleus of a disrupted 
galaxy (jAnnibali et al.|[2011al lbh: such observations indicate 
that the NGC 4449 may also have had more ancient ac- 
cretion activity, which may partially explain the significant 
amount of mass in the HI streams surrounding the system. 
Signatures of earlier a ccretion events in dwarf galaxies has 
also been presented by Geha et al. 1(|2005'), who find evidence 
for a counterrotating core in the elliptical dwarf galaxy NGC 
770 that they attribute to a minor merger event. 

Such observations clearly illustrate that LMC mass ob- 
jects do cannibalize smaller companions; however, the hall- 
marks of these encounters, such as faint tidal s treams, are 
challenging to observe. The IWilcots PrescottI ^004) con- 
clusions may thus indicate that the perturbing compan- 
ion has already been cannibalized, causing most Magel- 
lanic Irregulars to appear as isolated objects. Moreover, the 
Mbaryon/Mtotal ratio is a steep function of mass for these 
low mass systems; many dwarfs in the local group have ex- 
tremely large mass-to- light ratios. As such, a 1:10 total mass 
ratio companion may have a very discrepant stellar mass ra- 
tio, making the identification of such a companion challeng- 
ing observationally. 

In the context of the work presented here, Magellanic 
Irregulars are therefore key targets for deep HI and optical 
follow up observations as they are expected to be associated 
with tidal HI and stellar streams. Particular attention should 
be paid to Magellanic Irregulars with high current star for- 
mation rates, such as NGC 4449 and the LMC, which may 
point to ongoing tidal interactions with a low mass compan- 
ion. Furthermore, there is clearly need for future observa- 
tional and theoretical studies to better statistically quantify 
the frequency of interactions between LMC mass galaxies 
and smaller companions in order to assess the ubiquity of 
the theory presented here for the nature of Magellanic Ir- 
regulars. 

6.6 Assessment of the Models 

In this study we have explored the consequent large scale 
structure of the Magellanic System, the internal structure, 
kinematics of the LMC and the recent star formation histo- 
ries of the MCs in a first infall scenario; i.e. without strong 
tidal torques from the MW. We focus on two different mod- 
els for the interaction history of the LMC-SMC, one of which 
invokes a direct recent collision (Model 2), whereas in the 
other, the MCs never get closer than 20 kpc (Model 1). The 
ability of the presented models to reproduce key observed 
features of the Magellanic System are summarized in Table 

m 

Both models are able to reproduce the global large scale 
structure of the Magellanic System. Overall, however. Model 
1 provides better agreement with the properties of the Mag- 
ellanic Stream, whereas Model 2 provides significantly better 
agreement with the structure and kinematics of the LMC. 

While neither model reproduces every one of the fea- 
tures listed in Table [5] (and indeed the real answer is prob- 
ably somewhere in between the two presented models) it is 
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Table 5. Observed Properties of the Magellanic System and How the Models Fare 



Object 


Property 


Model 1 


Model 2 


Alternative 


Large Scale Structure 


A Leading Arm 


Yes 


Yes 






Location of Leading Arm 


No 


No 


Ram pressure 




A 150 degree Stream 


Yes 


Yes 






Stream location offset from orbit 


Yes 


roughly 


LC. : SMC disk inchnation 




Stream velocity gradient 


Yes 


roughly 


LC. : SMC disk inchnation 




Stream bifurcation 


roughly 


No 


Hydro instabilities 




Stream column gradient 


No 


No 


xxaiii jTieooLiie 










T QA/rr^ n-QC mace 




otream total mass 


No 


No 


Ram Pressure 5 
















Yes 


Yes 






in Rridp-p 


No 


Yes 






No ^F in ^trpam 


Yes 


Yes 




LMC 


Rotation Curve 


Yes 


Yes 






liflrQP'l" crQc nnri Q'i'pllpiT' 
v^iioct ^cxo aLL\jL otciicxi 


No 


Yes 






kinematics 










Offset bar 


No 


Yes 






Bar not seen in gas 


No 


Yes 






Warped stellar disk 


Yes 


Yes 






Elliptical stellar disk 


roughly 


Yes 






Current SFR 


No 


No 


SF details, feedback 


SMC 


Increasing SFR < Gyr 


No 


Yes 






Current SFR 


No 


No 


SF details, feedback 



Note. — SF stands for star formation. LC. stand for initial conditions. The column marked Alternative indicates other possible factors 
that may help fix discrepancies between the observations and the models. 



still rather remarkable that a single self- consistent model 
(namely Model 2) can simultaneously reproduce a large 
number of these features. Generally where the models fail 
(e.g. the column density gradient or location of the Lead- 
ing Arm) the likely missing ingredients are ram pressure 
stripping owing to the passage of the galaxies through the 
ambient hot gaseous halo of the MW or a detailed search of 
model parameters (e.g. initial gas mass or orbital parame- 
ters). Inconsistencies with the current SFR and recent SFHs 
are almost certainly a result of the star formation prescrip- 
tions employed and the lack of stellar feedback. 

A number of the discrepancies with the orbital model 
(e.g the SMC's position and velocity) can likely be addressed 
by a complete parameter search of, e.g., plausible mass ra- 
tios between the two galaxies (we chose a fixed mass ratio of 
1:10) and different orientations of the SMC disk relative to 
the LMC-SMC binary orbital plane. The inclination of the 
SMC's disk can dramatically alter the location and proper- 
ties of the Stream, e.g. a retrograde coplanar configuration 
would inhibit the formation of a stream entirely. It can also 
change the way LMC torques affect the SMC's motion. This 
is likely the explanation for why the Model 2 results do not 
reproduce the exact velocity field or location of the Stream. 

The Stream is also known t o be both spat i ally a nd kine- 
matically bifurcated, leading Putman et al.l (|2003h to de- 
scribe it as a "twisted helix". The Stream in our Model 
1 is not of constant column density along its width; be- 
tween Magellanic Longitudes of -50 and -100 the simulated 
stream appears to split into two high column density fila- 



ments. This bifurcation is largely a result of the rotation of 
the SMC's initial disk. Note also that because the Stream 
is seen in projection, the appearance of this bifurcation is 
highly dependent on the viewing orientation (which is model 
dependent). This is likely why similar structures are not 
seen as clearly in Model 2. But it should be noted that in 
the Gadget simulations the ISM is smoothed by an effec- 
tive equation of state, whereas in reality inhomogeneities 
should be present in the real ISM (e.g. molecular clouds). 
Stripping from a clumpier SMC ISM would lead to differ- 
ences in the distribution of the stripped material. We fur- 
ther expect that hydrodynamic instabilities, which are not 
modeled here, would capitalize on any initial density inho- 
mogeneities and augment this bifurc ation. It is clear from 
the significant observed turbulence (Nigra et al. 2010) and 
the exist ence of head-tail structure s in cloudlets within the 
Stream (|Putman et al.|[2003lj201lb , that hydrodynamic in- 
stabilities are shaping its internal structure. Gas drag will 
operate differently on these clumps depending on their den- 
sities, naturally leading to a velocity bifurcation as well. 

Although strong perturbations in the LMC disk are ex- 
pected from a close encounter with a 1:10 mass ratio com- 
panion, without a large parameter study of impact param- 
eters for the collision and mass ratios, it is unclear how 
generically these results can be applied to Magellanic Ir- 
regulars. A future study of these parameters will assess the 
robustness and longevity of the presented asymmetric struc- 
tures (i.e. offset bars and one-armed spirals). In particular, 
it remains unclear as to whether a direct collision is always 
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required to produce such structures. Regardless, the exis- 
tence of Magellanic Irregulars with close pairs connected by 
gaseous bridges is certainly suggestive of a link to a similar 
interaction history as that of the MCs. 

It is worth pointing out that while it has been spec- 
ulated that many of the features listed in Table [5] were 
directly related to interactions with the SMC, most of these 
links have never been illustrated by numerical simulations 
which self-consistently reproduce the global large scale fea- 
tures of the system. In this light, particular successes of the 
simulations presented in this work include the reproduction 
of a warped, off-centered bar that is neither detectable in 
the gaseous disk nor actively forming stars while simultane- 
ously forming a Bridge, Leading Arm and trailing 150 degree 
long Stream. Furthermore, we have not yet discussed the re- 
markable agreement of the simulated SMC model with its 
observed structure and kinematics from this same Model 2. 
These results will be outlined in a subsequent paper. 



7 CONCLUSIONS 

We have explored two models for the possible interaction 
history of the LMC and SMC in an effort to simultaneously 
reproduce both the large scale gaseous distribution of the 
Magellanic System and the internal structure and morphol- 
ogy of the LMC. Here we summarize our findings for the 
Magellanic System and the implications for the study of 
Magellanic Irregular galaxies more generally. 

7.1 Conclusions for the Magellanic System 

The resulting kinematics and structure of the LMC strongly 
favor a scenario in which the MCs have recently (100-300 
Myr) experienced a direct collision (our Model 2). Orbital 
models where the SMC never gets closer than 20 kpc to the 
LMC (e.g. Model 1) are able to reproduce the large scale 
structure of the Magellanic System, but poorly match the 
LMC's internal properties. In particular, without a direct 
collision with the SMC, the LMC would be better described 
as a normal, symmetric barred spiral galaxy. An upcoming 
paper will also illustrate that the observed internal kinemat- 
ics and morphology of the SMC are also better described by 
this same collision model (Besla et al., in prep). 

This study illustrates that, surprisingly, the LMC's disk 
can maintain a fairly smooth stellar velocity field despite a 
direct collision with the SMC and that such a scenario can 
explain a number of observed features of the LMC disk: 

(i) The gas and stellar kinematic centers of the LMC disk 
are coincident; however they are offset from the photometric 
center (Cole et al. 2005). 

(ii) The old stellar disk of the LMC is thick and warped. 
The warping at the e dges may give it a flared appearance 
(jAlves Nelsonll2000h . 

(iii) The bar is warped relative to the LMC disk plane by 
ten degrees and is off-center relative to the dynamical center 
of the gaseous disk. While an offset bar in the LMC has been 
suggested by numerous aut hors as being the result of col- 
lisions with the SM C (e- g -, ISubramaniam SubramanianI 
I2OO9I : ISubramaniam" 2003; Bekki k Chiba 2007') this is the 
first time it has been modeled self-consistently. 



(iv ) The bar is not seen in the gas distribution (|Kim et al.l 
llQQSh or as a site of on- going star formation, likely because 
the bar is warped out of the plane, inhibiting efficient gas 
funneling. 

(v) Gas eous "arms" similar to th ose seen in HI maps of 
the LMC (iNidever et al.ll2010l . l2008h are stripped out of the 
LMC by the SMC in the direction of the Magellanic Bridge. 
This "arm" was formed during a violent collision and may 
have signatures in polarization maps of the LMC's magnetic 
field (Mao et al. in prep). This result also implies that there 
should be a metallicity gradient along the Bridge, increasing 
towards the LMC owing to contamination by LMC gas. 

(vi) A one-armed spiral is induced in the LMC's disk and 
is a site of on-going star formation. 

(vii) Stellar debris from the SMC is expected in the same 
field of view as the LMC disk. These stars will have differing 
kinematics signatu res from the local LMC disk velocity field 
(jOlsen et ahlbOllh and may be the source of microlensing 
events towards the LMC (Besla et al. in prep). Such tidally 
stripped stars are not expected in a pure ram pressure strip- 
ping model for the Stream. 

(viii) The gaseous Bridge that connects the two galaxies is 
the site of on-going star formation, where the SFR increases 
along its length towards the SMC. 

(ix) The Leading Arm extends >70 degrees ahead of the 
MCs and the Magellanic Stream extends 150 degrees behind 
them, reproducing the full extent of the Magellanic System. 
The Leading Arm is a younger structure than the Stream 
since it formed as a tidal tail during the SMC's most recent 
orbit about the LMC. 

These listed properties are a direct consequence of 
interactions with the SMC^ con firming the suspicions of 
Ide Vaucouleurs FreemanI (|l972h that the LMC's peculiar 
morphology does not owe to interactions with the MW. 

7.2 Implications for Magellanic Irregulars 

There is potentially much to learn about Magellanic Irreg- 
ular galaxies as a class by studying the Magellanic System 
in detail. In this work we have shown that the off-center bar 
and one-armed spiral arm of the LMC may be a product 
of close encounters with its smaller companion, the SMC. 
We conclude that interactions between a massive dwarf and 
a smaller companion, even one 10 times smaller in mass, 
can significantly alter the morphology of an otherwise nor- 
mal looking, low-mass spiral galaxy. This study thus indi- 
cates that dwarf-dwarf galaxy interactions can be important 
drivers of their morphological evolution, without relying on 
interactions with a massive host. 

Given that off-center bars and one-armed spi- 
rals are common characteristics of Magellanic Irregulars 
(jde Vaucouleurslll964f ), it is possible that such galaxies cur- 
rently have or once had small companions. Known examples 
of Magellanic Irregulars with a close companion are thus 
prime candidates for follow up HI surveys to map the ex- 
tended HI distribution, as the bridges that connect them 
likely have tidal tail counterparts. We argue here that if 
such systems were accreted by a massive spiral, they could 
form an analogous Magellanic System, i.e. two interacting 
dwarfs surrounded by an extended HI complex. A poten- 
tial example could be the interacting pair of MC analogs. 
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NGC 4485/4490, which are surrounded by an extensive 
HI envelope but are not in close proximity to a massive 
host (Cle mens et al.|[l998l ). The recently identified stellar 
stream about the Magellanic Irregular galaxy NGC 4449 
(|Martinez-Delgado et alJEoilh gives further credence to this 
theory. 

Quantifying the efficiency with which such dwarf-dwarf 
tidal interactions can remove gaseous material is directly 
relevant to questions about how dwarf galaxies lose their 
gas and how gas is supplied to more massive galaxies. The 
frequency of such dwarf-dwarf encounters has not yet been 
quantified observationally or theoretically. As such, the rel- 
ative importance of such encounters remains to be deter- 
mined. However, by reproducing the large scale structure 
of the Magellanic System in a LMC-SMC tidal scenario, we 
have shown here that at least 50% of the original gas budget 
of a small dwarf can be easily removed by tidal interactions 
with a larger dwarf companion, making such encounters a 
potentially important mode of baryon loss on these mass 
scales. 

While the accretion of an LMC+SMC analog by 

a MW type host is a relatively rare event today 

(|Bovlan-Kolchin et ahlboill : iLiu et al.ll201lh . modern mod- 
els for the cosmological assembly of galactic- scale structures 
suggest that the MW halo formed from the earlier accre- 
tion and disruption of LMC-mass objects (| Stewart et al.l 
I2OO8I ). As such, simulations of isolated Magellanic Irregu- 
lars (LMC analogs) with low mass companions and the di- 
rect comparisons of such simulations to o bserved analogs, 
such as NGC 4027 (|PlK)okun et al.l ll992') and NGC 3664 
(|Wilcots Prescottl l2004l), provides a logical testing ground 
for our understanding of the morphological, kinematic and 
chemical evolution of the fundamental building block of a 
MW-type galaxy. 
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APPENDIX A: THE ROLE OF RAM 
PRESSURE STRIPPING IN THE FORMATION 
OF THE STREAM 

Evidence for the existence of a hot gaseous halo compo- 
nent o f our MW comes from absorption lines towards distant 
AGN (|Sembach et al.ll2003h . The MCs are therefore known 
to be moving through some ambient medium, and so the ef- 
fects of ram pressure stripping and gas drag likely play some 
role in explaining the gas distribution of the Magellanic Sys- 
tem. Our contention is that the impact of ram pressure is to 
modify the properties (velocity, location, mass) of the gas 
distribution, but does not in itself explain the existence of 
the Stream, Bridge and Leading Arm Feature. 



To model the impact of ram pressure from a physical 
ambient medium representing the MW's gaseous halo with 
either an SPH or Eulerian grid-based code would require 
extremely high numerical resolution, as the gaseous stream 
is very diffuse. Instead, we have created a toy simulation, 
where the Gadget-3 code was modified to include a ram pres- 
sure acceleration term that is applied uniformly to all gas 
particles and operates opposite to th e direction of n iotion 
of those gas particles. Following the ICunn GottI (|l972l l 
prescription, the ram pressure experienced by a gas particle 
moving through the hot gaseous halo of the MW (density, 
Phot), at a speed Vg^i is expressed as follows: 



^ PhotVg 



(Al) 



The acceleration experienced by a gas cloud of surface 
density Sgas owing to ram pressure is: 



— -fram/5^g 



(A2) 

Following IVollmer et al. I (|200lh . we can describe an 
HI gas cloud as having a characteristic co lumn density of 
Nm = 7.5 X 10^°cm-2 (ISanders et al.|[l985l ). Thus, Sgas - 
Nhi X TTiHi- For a given velocity component Vi=x,y,z, the 
corresponding acceleration applied in that direction per 
timestep is: 



ai 



^hot 2Vi_ 

"Nhi'' v' 



(A3) 



where the volume density (phot) has been replaced by 
the number density of the MW's ambient halo medium 
(rihot) multiplied by the mass of the average particle in 
the halo, which we have assumed is ^ mm- The quantity 
rihot is largely unconstrained , although an upper limit of 
5 X 10~^ cm~^ is estimated ( IRasmussen PedersenllioOll : 
iMalonev k Bland-Hawthornlll999h We consider three val- 
ues for this parameter: 5 x 10~^, 10~^ and 5 x 10~^ cm~^. 
We further assume that this density is constant all the way 
to R200 = 220 kpc of the MW. The ram pressure accelera- 
tion is applied to all gas particles as soon as the MCs enter 
within R200 of the MW. Note that this acceleration will not 
be constant as a function of time since the velocity of the 
particles changes as the MCs fall towards the MW. 

We stress that this is a crude model, as there is no ac- 
tual ambient gas density present and the force is applied uni- 
formly to all particles. In reality, gas in the inner parts of the 
disk should be shielded by those on the outside. Moreover, 
other hydrodynamical instabilities, such as Kelvin-Helmoltz 
and Ray leigh- Taylor, will not be modeled in this set-up and 
we will therefore underestimate the amount of gas loss the 
galaxies may incur. This toy model will, however, give us 
a rough idea of how the position of the simulated Leading 
Arm, Bridge and Stream will evolve over time. 

In Figure lAll and IA2I the gas column density of the 
Magellanic System is mapped in Magellanic coordinates for 
Model 1 and Model 2, respectively. As the background den- 
sity is increased, the location of the Stream cha n ges, a nd 
the gas begins to trace the orbit. In iBesla et"al1 (|2007f ) it 
was shown that the location of the LMC past orbit deviates 
from the location of the Stream on the plane of the sky. This 
result was further shown to be insensitive to the MW model 
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Figure Al. The gas column density in Magellanic Coordinates along the Magellanic System for Model 1 with a toy model ram pressure 
acceleration. The ambient hot halo gas density increases from top to bottom, as marked. Ram pressure causes the simulated stream to 
trace the past orbits of the Clouds on the plane of the sky; this is opposite to the location of the true Stream (indicated by the white 
line). The column density along the simulated stream also changes as the ram pressure becomes more important. 



and robust within 3a of the measured proper motions. The 
fact that ram pressure stripping works to align the stripped 
material and the past orbits suggests that it is not the main 
formation mechanism for the Magellanic Stream. 

The Leading Arm also changes location and structure 
significantly in the toy model, particularly in Model 1. Fac- 
ing gas densities larger than 10 ~^ cm~^ it disappears en- 
tirely. Although the Bridge disappears in Model 1 at densi- 
ties larger than 5 x 10~^ cm~^, it remains a strong feature in 
Model 2 at all densities. This is likely because, in Model 2, 
the Bridge is a young feature that formed in the most recent 
collision between the MCs. It has thus not had enough time 
to experience significant ram pressure effects. In Model 2 the 
leading arm also gets closer to the MW, falling to a line-of- 
sight distance of 30 kpc. These models do not predict that 



components of the Leading Arm should be interacting with 
the gaseous disk of the MW, although the re are claims that 
such a situation has been observed (McClure- Griffiths et al.l 
I2OO8I I. 

Note that as the ram pressure increases, the SMC's orbit 
changes; the SMC is clearly not in the same location in the 
lower panel of Figures lAll and IA2l as in the respective upper 
panel. This occurs because ram pressure decreases the veloc- 
ity of the SMC. The LMC, on the other hand, is too massive 
for its motion to be affected. A strong ram pressure head- 
wind is thus likely at odds with the high relative velocity 
observed between the MCs (^^100 km/s, .Kallivavalil et al.l 
l2006bh . 

The structure of the LMC gas disk changes even with a 
mild ram pressure headwind. The gas disk is rotating clock- 
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Figure A2. Same as Figure [All but for Model 2. 
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wise, and so the lower half is rotating into the head wind 
and gets stalled. Gas therefore builds up in the lower left 
corner of the LMC's disk. This is in f act where the 30 Do - 
radus star-forming region is situated. I de Boer et al. I([l99i) 
suggested that the HI overdensity seen in the South-East is 
a direct result of the interaction between the rotating LMC 
gaseous disk and a headwind from the ambient medium. 
Here, we can indeed see that this process occurs generically 
in both Models 1 and 2. The LMC's gaseous disk is observed 
to be truncated (i.e. it is not as extended as the stellar disk). 
The simulated LMC gas disk is truncated if the ambient gas 
density is at least of order 5 x 10 ~^ cm~^. 



Note that the column density along the Stream also 
changes as ram pressure become more important. The col- 
umn density increases in the regions closest to the MCs, 
indicating that ram pressure may be the solution for the 
mismatch between the observed maximal column density 



gradient along the Stream and the simulation results pre- 
sented in Figure [SI 

From this simple toy model we can conclude that ram 
pressure and other hydrodynamic instabilities will change 
the mass budget in the various components of the system 
and change their locations on the plane of the sky. In partic- 
ular, ram pressure works to align the Stream with the past 
orbits, contrary t o expectations from the proper motions 
(jBesla et al.ll2007l l . We estimate that if the background halo 
gas densities are in excess of 10~^ /cm^ the gas distribution 
of the simulated Magellanic System will be irreconcilable 
with observations. This estimate is in accord wit h observa- 
tional upper limits (jRasmussen &; PedersenllioOll ). 

A recent study bv lDiaz (|2011bl l suggests that 

ram pressure at the am bient densities required by the 
iMastropietro et"aL 1 (|2005l ) study (5 x 10"^ /cm^) would be 
too large for the survival of the Magellanic Stream. In a first 
infall scenario, where the Magellanic System would be inter- 
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acting with the ambient halo medium for a shorter period of 
time (< 1 Gyr), we find that the Stream can indeed survive 
such densities and that ram pressure stripping hkely plays a 
very important role in shaping the Magellanic Stream (par- 
ticularly the Leading Arm) and increasing its mass budget. 
It should also be noted that Faraday rotation has been de- 
tected in at least one high velocity cloud in the Leading Arm 
(|McClure- Griffiths et al. 2010); magnetic fields may protect 
the cloudlets from evaporation and hydrodynamic instabili- 
ties, increasing the chance of survivability of the Stream in 
the face of a ram pressure headwind. 



APPENDIX B: THE ROLE OF STELLAR 
FEEDBACK IN THE FORMATION OF THE 
STREAM 

Ther e is evidence of stefi ar feedback in the vicinity of the 
MCs. lLehner et aP (| 20091 ) have detected high velocity clouds 
(HVCs) between the MCs and the MW, moving at a velocity 
with respect to the local standard of rest (visr) as high as 150 
km/s. Using FUSE, they find on average that these HVCs 
have metallicities of [OI/HI] = -0.51 ie masses of 
(0.5-1) xlO^M©, although these HVCs are predominantly 
ionized. Such material is thus, on averag e, enriched rela- 
tive t o the LMC and SMC. Furthermore, Lehner &: Howkl 
(|2007l l have detected a highly ionized (OVI, CIV, SilV, NV) 
corona of high velocity gas surrounding th e LMC, suggestive 
of outfiows. Stavelev-Smit h et al.1 (l2003bl ) also report simi- 
lar observations, and find some of this high velocity gas to 
be projected on HI voids in the LMC. Data from FUSE in- 
dicates that the SM C appears to also have an OVI corona 
(|HooDes et al.ll2002l ). These detected structures are in close 
proximity to the MCs, and so may be indicative of mass 
loss processes that are currently ongoing, rather than a long 
duration process that would be required to build the Mag- 
ellanic Stream^ 

Recently, iNidever et al. (2008) found a coherent visr 
gradient from the LMC along one of the filaments in the 
Stream and in the Leading Arm. They also claim to detect 
sinusoidal velocity patterns in the Stream, which they in- 
terpret as signatures of the LMC's disk rotation. They thus 
conclude that as much as half of the mass within the Stream 
originates from the LMC and that the sinusoidal velocity 
pattern indicates that the bulk of this material is emanat- 
ing from the south-east HI overdensity (SEHO) region of 
the LMC's HI gaseous disk as a stellar outfiow over the past 
1.7 Gyr. The idea that the Magellanic Stream was formed 
from stellar outfiows was also suggested by aOlang (,20o3 ) . 
The LMC is pock-marked with giant superbubbles, indicat- 
ing locations of strong stellar feedback/ winds. If the column 
density were constant across the disk, then each superbub- 
ble would have originally contain e d roughly 10^ Mc^) w orth 
of material (|Kim et al.ll 19981 . 1 19991 ). iNidever et HI (|2008l ) de- 
termine that 2-3 superbubbles losing 10^ Mq worth of mate- 
rial every 10 Myr over 1.7 Gyr would be sufiicient to explain 
the mass budget of the Stream. 

A stellar feedback model for the Stream is attractive in 
that it does not rely on strong MW tides to remove mate- 
rial from the LMC, making it consistent with a first infall 
scenario (although the leading component would still be an 
issue) . 



However, it is unlikely that outfiows generated by SNe 
feedback would be sufficiently energetic to be unbound from 
the LMC's gravitational potential. In our first infall LMC 
models {Mlmc = 1.8 x 1 0^^ Mp) , the e scape speed is of or- 
der 250 km/s at 10 kpc. iMartinI (|2005l ) finds that terminal 
velocities expected from galaxies with maximum circular ve- 
locities of 100 km/s should also be of order 100 km/s, well 
below the escape speed at 10 kpc. Moreover, galaxies with 
star formation rates as low a s 0.1M(;)/yr (i.e. the va lue in 
the 30 Doradus region todav. [Harris Zaritskvl [20091 ) gen- 
erally have outfiow velocities less than 30 km/s (figure 6 of 
Martin 2005). 

However, it is certainly possible that a ram pressure 
headwind can exploit the stellar feedback processes to aid 
in the removal of material from the deepest parts of the 
LMC's potential. This theory is testable by observations of 
the metallicity in the Stream: not only do the LMC and SMC 
have different metallicities, but any feed back scenario will re- 
sult i n the removal of enriched material ([Mac Low Ferraral 
1 19991 ). 

Metallicities for various components of the Magel- 
lanic System are summarized in Table IB1[ The cur- 
rent day LMC's m etalhcity is [0/H]lmc = -0.34 zb 0.06 
(iRussell Dopital Q,992 , updated to the latest solar abun- 
dances by Fox et al. 2010). B ut chemical enrichment mod- 
els of IPagel Tautvaisiend ([l998) suggest that 2 Gyr 
ago the LMC's metallicity could have been as low as 
[0/H]lmc ^ -0.5 (their Figure 3). 1-2 Gyr i s the relevant 
timescale for the formation of t he Stream (jNidever et al.l 
'2008'; 'Gardiner k Noguchi' "l996); gi ven the measured Ha 
emission along its length (^Weiner Williamslll996 ) and the 
corresponding expected abla tion timescale for the neutra l 
HI component of the Stream ([Bland-Hawthorn et al.ir2007l ), 
it is unlikely that the Stream could have survived for much 
longer. 

Recently, I Fox et aP (|2010l ) have determined the oxygen 
abundance for a region near the tip of the Stream from ab- 
sorption lines towards the background quasar NGC 7469. 
Observed [OI/HI] abundances are a close indication of the 
true oxygen abundance, [0/H], s ince oxygen is not s trongly 
depleted onto interstellar dust (| Jensen et al.l l2005l ) . They 
find: 
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HI 



1 /' N(oi) A 



log 



log 



10 



14.32±0.04 



]^Q18.63±0.13 

1.00 ± 0.05(stat) ± 0.08(syst) 



(Bl) 



3.31) (B2) 
(B3) 



In an outfiow or ram pressure stripping scenario the 
metallicity of the Stream must be at least that of the origi- 
nal ISM. These measurements make it improbable that the 
Stream could have solely originated in the LMC, regardless 
of the formation mechanism. The SMC interstellar oxygen 
abundance is a better match to these observations (see Table 
IBJ,D, especially since it was also less enriched 1-2 Gyr ago 
(Pagel k Tautvaisiend Il998l ). Interestingly, such low abun- 
dances are also measured in th e Bridge connecting the MCs 
f Table [ml ILehner et al.ll2008l ). suggestive of a common ori- 
gin for both structures. 

Closer to the MCs, the metallicity measurements de- 
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rived from absorption lin es towards background quasars are 
larger by a factor of 2-4 (|Lu et al.l ll998': Gibson et al.l l2QQ0l : 
ISe mba ch et alj I2OOII) . But these other measurements also 
have much larger error bars, as they use tracers that are 
more sensitive to ionization corrections than the OI/HI ra- 
tio. If we take an upper limit of a factor of 4 incr ease in the 
oxygen abundance in the Stream relative to the IFox et ahl 
(I2OIO) values, we can assess whether this material could have 
once (^-^2 Gyr ago) originated in a s tellar wind from th e 
LMC. We follow the methodology of iMartin eraP (|2002h . 
who estimate the expected metallicity of outflows from the 
starbursting dwarf NGC 1569. The total mass of the ejected 
wind Mw is given by 



(B4) 



where Mej is the mass in SNe ejecta and x is the mass load- 
ing of the wind. We consider first the required mass loading 
from the ISM of the LMC needed to dilute the metallicity 
so as to not violate the upper limits for the Stream. The 
oxygen abundance of the wind can be expressed as. 



Mej(Zo,sN + xZo,ism) 
Mw 



(B5) 



Written in terms of the solar values, the abundance of 
the ISM 2 Gyr ago is Zo,ism = 10"°-^ = 0.32. Following 
iMartin et all (|20Q2l ). the IMF-averaged metaflicity of SNe 
ejecta is Zo,sn ~ 8 times solar for oxygen. We choose three 
values for Zo,sn = 4, 8 and 12. We can rewrite equation IB5I 
to solve for x ^ mass independent way using equation lB4l 



(B6) 



Given that we know the current upper limit for the 
metallicity of the wind, Zo,w = 10~°"^ = 0.4, then the re- 
quired mass loading is x(Zo,sn = 4,8, 12) = (45,95, 145). 

Mass loading in excess of factors of 10 are not observed. 
The wi nd from M82 is expec ted t o be mass loaded by a factor 
of 3-6 (I Suchkov et allll996l ) and IMartin et al.1 (|2002l ) find a 
value of X = 9 to be fav ored for NGC 1569. Furthermore, 
numerical simulations bv lMac Low &: Ferraral (| 19991 ) find it 
diflficult to accelerate any ambient cool gas that is swept up 
with the expanding outflow; instead, much of that material 
remains bound to the galaxy. 

Ram pressure stripping might be the key missing ingre- 
dient that could help entrain ambient gas and explain such 
high mass loading factors. Assuming the addition of ram 
pressure stripping results in mass loading in excess of a fac- 
tor of 45, then the resulting total wind mass (equation IB4p 
over a timescale of 1.7 Gyr (the lifetime of the Stream ac- 
cording to the Nidever et al. 2008 model) would be of order 
While this appears to be more than enough 
material to explain the total amount of neutral HI observed 
in the Stream (5 x 10^ M©; Table[4]), it would imply that the 
entire Stream should be metal enriched, contrary to obser- 
vations. This scenario thus over-predicts the mass budget of 



2 Mej = FM^Tms is determined assuming a star formation rate 
of M* =0.1 M0/yr (the current rate in 30 Doradus), a S Ne forma- 
tion rate of F =0.12 determined from StarburstQQ (|Leitherer et al.l 
[19991 with a Kroupa IMF, over a timescale of Tms =1-7 Gyr 



the Stream and requires mass loading factors that are more 
extreme than observed. 
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Table Bl. Oxygen Abundances 



Object 



[O/Hl^ 



Reference 



LMC (today) 
LMC (2 Gyr ; 
SMC (today) 
Stream 

Stream (Upper Limit^ 
Bridge 

HVCs near LMC 



-0.34 ± 0.6 

-0.5 
- 0.66 ±0.1 
-1.0 ±0.13 

-0.4 
-0.96ig:}? 



Russell &: Dopita (1992); Fox et al. (2010 



Pagel Tautvaisiene (1998) 
Russell Dopita (1992) : Fox et al. (2010) 
Fox et al. (2010) 
Gibson et al. (2000) : Lu et al. (1998); Se mbach et al. (2001) 
Lehner et al. (1008 ) 
Lehner et al. (2009) 



^ Log of the oxygen abundance relative to solar, as defined in equation IB2I 

^ Upper limit is defined by increasing the observed oxygen column density of the Stream by a factor of four 
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